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Assuming the validity of random matrices for describing the statistics of a closed chaotic quantum 
system, we study analytically some statistical properties of the S-matrix characterizing scattering 
in its open counterpart. In the first part of the paper we attempt to expose systematically ideas 
underlying the so-called stochastic (Heidelberg) approach to chaotic quantum scattering. Then 
we concentrate on systems with broken time-reversal invariance coupled to continua via M open 
channels; a — 1,2, ..,M. A physical realization of this case corresponds to the chaotic scattering 
in ballistic microstructures pierced by a strong enough magnetic flux. By using the supersymme- 
try method we derive an explicit expression for the density of S-matrix poles (resonances) in the 
complex energy plane. When all scattering channels are considered to be equivalent our expression 
describes a crossover from the distribution of resonance widths (regime of isolated resonances) to 
a broad power-like distribution typical for the regime of overlapping resonances. The first moment 
is found to reproduce exactly the Moldauer-Simonius relation between the mean resonance width 
and the transmission coefficient. Under the same assumptions we derive an explicit expression for 
the parametric correlation function of densities of eigenphases 9a of the S-matrix (taken modulo 
2-k). We use it to find the distribution of derivatives Ta = d9a/dE of these eigenphases with respect 
to the energy ("partial delay times" ) as well as with respect to an arbitrary external parameter. 
We also find the parametric correlations of the Wigner-Smith time delay Tw{E) — -jj dOg/dE 
at two different energies E — Q/2 and E + i}/2 as well as at two different values of the external 
parameter. The relation between our results and those following from the semiclassical approach as 
well as the relevance to experiments are briefly discussed. 
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I. INTRODUCTION. 



Chaotic scattering has been a subject of a rather intensive research activity during the last decade, both theoreticaUy 
(see reviews and experimentahy ^|-^2j.This phenomenon is encountered in a variety of physical systems ranging 



from atomic nuclei P, [L3||l4| , atoms |l5|-17|] and molecules 19 1 to mesoscopic ballistic devices P,p()[ and microwave 
cavities . The most fundamental object characterizing the process of quantum scattering is the unitary 5*— matrix 

relating the amplitudes of waves incoming onto the system and the amplitudes of scattered (outgoing) waves. Because 
of the chaotic nature of the underlying scattering dynamics the 5— matrix characteristics behave in an irregular way 
when parameters of either incoming waves (e.g. energy) or of the scattering region (e.g. the form or strength of the 
scattering potential, the strength of the magnetic field through the ballistic microstructure, etc.) are slightly changed. 
Because of this fact it seems to be most adequate to describe such a behavior in terms of some statistical measures: 
distributions and correlation functions. 

At present, there are two complementary theoretical tools employed to calculate statistical characteristics of open 
quantum systems exhibiting the phenomenon of chaotic scattering. These are the semiclassical [plj^,^, ^,^ and 
the stochastic approaches the relation between both methods being in some detail discussed in |25|. The 

semiclassical approach operates with the genuine microscopic Hamiltonians and allows for treating particular systems 
with full account of their specific features. The starting point for such an approach is a representation of ^-matrix 
elements in terms of a sum over the classical periodic orbits, the method going back to works by Gutzwiller |26[ | and 
Balian and Bloch |^ . The statistical characteristics are sampled usually over some range of energies or changing the 
system parameters. 

It is however known, that the majority of (both closed and open) quantum chaotic systems of quite different 
microscopic nature shows a great degree of universality in their properties on the appropriate energy scale. More 
precisely, the statistical characteristics of closed systems turn out to be independent of the microscopic details when 
sampled on the energy intervals SE large in comparison with a mean separation between two adjacent levels A, but 
smaller then the energy scale Ec — h/te, with standing for the relaxation time necessary for the classically chaotic 
chaotic system to cover uniformly the constant energy shell [|2^ . Because of this universality one achieves the correct 
description of the properties of such systems [ p9| by exploiting the similarity with ensembles of large Gaussian random 
matrices H of the size N x N characterized by the following probability distribution: 

Vpocexp-^TrH^ (1) 

where the matrices are considered to be real symmetric {13 — 1, Gaussian Orthogonal Ensemble: GOE), Hermitian (/3 — 
2, Gaussian Unitary Ensemble: GUE) or consisting of real quaternions {(3 = 4, Gaussian Symplectic Ensemble: GSE). 
Ensembles with P = 1 {(] — 2) serve to describe spectra of closed quantum chaotic systems systems with preserved 
(broken, e.g by applied magnetic field or by Aharonov-Bohm magnetic flux) time-reversal invariance, correspondingly. 
At last, the ensemble corresponding to /3 = 4 describes systems with preserved time-reversal invariance displaying 
strong spin-orbit scattering which should be taken into account. 

Properties of all these ensembles were studied long ago [ ^Opl| ]. The mean level density Vsc{E) ~ {^{E)) in the limit 
iV — !■ cx) is given by the so-called Wigner semicircle law: 

Vsc{E) ^ (^^Tr6{E - H)^ = ^V^^ (2) 

where the angular brackets stand for the ensemble averaging. The radius of the semicircle is equal (in chosen normal- 
ization) to Esc = 2, so that the average spacing between eigenvalues is 4/A^ while the local spacing around the point 
E is A(S) = {NuscT^. 

The mean level density is the simplest quantity characterizing the spectrum of any system. It is not very informative 
from a physical point of view since it is insensitive to the fine structure of the spectrum. It is also the same for all 
universality classes. Actually, no real physical system is known where the mean level density follows the semicircle 
law, Eq.(|2|). 

In contrast, the two-point spectral correlation function: 

R2{uj) ^ /^^{E){v{E ~VL)v{E + n)) ~l (3) 

where w = 27rr2/A(i?) is known to be universal, i.e. independent of the microscopic details and has the same form 
both for generic chaotic systems and for the Gaussian Ensembles of definite symmetry. The same universality extends 
to other spectral properties, such as the nearest neighbors spacing distribution, etc. [p9[. 
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Despite the apparent success in the exploitation of random matrix resuhs for describing spectra of quantum chaotic 
systems of different nature p9| , ^ it was a long standing problem to justify the validity of such an approach mi- 
croscopically. Some insight was achieved within the semiclassical approach long ago by Berry Very recently 
Muzykantsky and Khmelnitsky and especially Andreev et al. managed to find a way to prove this conjecture 
by a nontrivial combination of ficld-thcoretical and semiclassical ideas. In parallel, traditional semiclassical methods 
were also significantly improved |35| . These results put applications of random matrices for the description of universal 
features of closed chaotic systems on a firm ground. 

Provided the properties of a Hamiltonian 7ii„ for a closed chaotic system are specified, one can consider its open 
counterpart and work out the S— matrix by standard methods in the theory of quantum scattering p6|-^l[ . As the 
result, the scattering matrix is expressed in terms of both the Hamiltonian 7ii„ and matrix elements describing the 
coupling of the internal motion to "open channels" i.e. the states of the system asymptotically far from the chaotic 
region. Correspondingly, one can try to extract the statistics of S*— matrix inherited from the mentioned universal 
"random matrix" properties of Hin [ p2[ . 

In principle, it is far from being obvious that the coupling to continua does not wash out the universal features. 
The key observation (made long ago in the context of nuclear physics, see e.g. is that typically there are two 

well-separated time stages associated with the scattering process: an immediate "prompt" response (so-called direct 
processes) and a delayed, or equilibrated response associated with the formation of long-living states, or resonances. 
In the energy domain direct processes are described by smooth S*— matrix characteristics averaged over a large energy 
interval. Such characteristics must be, of course, highly non-universal and are determined mainly by system-specific 
boundary conditions on the boundary of the scattering region. At the same time, resonance response happening on 
much shorter energy scale manifests itself in a form of a random signal on top of the smooth averaged characteristics. 
Formation of the long-living resonances is intimately related to the internal dynamics inside the scattering region. It 
is natural to expect that the universal features of the chaotic quantum dynamics will be reflected in the universal 
statistical characteristics of the 5— matrix on the "resonance" energy scale, as long as the characteristic times(e.g. 
measured by inverse widths of the resonances) will be much longer than the time scale of the direct response. 

To find some adequate description of these universal features one can substitute the Hamiltonian Tim by the 
Gaussian random matrix, Eq.(0). This way was pioneered by Verbaarschot, Weidenmiiller and Zirnbauer p2[ who 
calculated the correlation function of S— matrix elements at two different energies for arbitrary number M of open 
channels satisfying M <^ N, with N being the total number of resonances. It was indeed found that the N x M 
matrix elements describing the coupling of the internal region to open channels enter the final result in the form of 
only M simple combinations, the so-called "transmission coefficients". In full agreement with the "two distinct time 
scale" picture, these coefficients just measure the portion of the flux in a given channel which is not reflected back 
immediately, but penetrates the interaction region and participates in the formation of the long-living resonances 
1^,^. The approach developed in (following ||2j] we will call it "Heidelberg approach" henceforth) turned out 
to be very fruitful and serves as a case study for all further development in the field. 

One can also try to make use of the expected universality directly on the level of S'-matrix without any reference 
to the system Hamiltonian. Such a method was developed in great detail in a series of papers by Mello and co- 
workers The probability density for the whole S'-matrix can be obtained if one makes the assumption of 
minimal information content of such a distribution respecting the requirements of S'-matrix unitarity, analyticity and 
constraints imposed by absence or presence of the time-reversal invariance. Provided all the system-specific relevant 
information is encoded in the value of the average S'-matrix (S) the joint probability P{S)dfj,{S) is given by: 

r ^ X 1 (/3M+2-/3)/2 

„ „ det(/-(S)(St)) 

I det(/-S(St)j \im+2-P) 

with the following measure d^ii^{S) : 

M 

dfif,{S) oc II I e'^» - e'^" I'' n dOadU (5) 

a<b a=l 

where 9a', a = 1, ...,M are eigenphases of the S— matrix, the volume element dU is generated by the corresponding 
eigenvectors, /3 = 1,2,4 as before and M is the dimension of the S— matrix equal to the number of open channels. 
The distribution P{S) is called the Poisson's kernel. 

For the particular case (S) = the distribution Eq.(^) just coincides with the measure d^f3{S). Such distributions 
were considered long ago by Dyson and known as the Dyson Circular Ensembles pO[ . They were found to describe 
very satisfactorily the S'-matrix statistics for some realistic models of the chaotic scattering [Esl . The general Poisson's 
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kernel, Eq.(Q) was verified as well |2^j2J| and proved to be a very useful tool to predict fiuctuations of transmissions 
through ballistic microstructures. It is natural to expect that the same distribution can be actually derived from the 
Heidelberg approach. It turned out that the problem is quite involved technically, however. In his insightful paper 
p6| Brouwer succeeded to derive the Poisson's kernel distribution assuming that the Hamiltonian 7ij„ is taken from a 
quite specific Lorentzian ensemble of random matrices. Since the spectral properties of the latter ensemble and those 
of Gaussian matrices, Eq.(|l|), are identical as long as the matrix size N oo, one expects that the equivalence of the 
two approaches can be shown for this generic case as well. 

If one wishes to study the dependence of the 5— matrix on external parameters without explicitly considering 
the system Hamiltonian, one should make some additional statistical assumptions beyond the minimum information 
approach. One possible way is to simulate such a dependence by a kind of " Brownian motion" in the corresponding 
S— matrix space p7| . It turns out, however, that the Brownian motion picture is in disagreement with the results 
obtained starting from the Heidelberg formalism. Therefore, the Heidelberg approach seems to be the only consistent 
stochastic method when we are interested in the parametric variations of the S"— matrix characteristics. An example of 
such kind of calculation can be found in |Q . Another important advantage of the Heidelberg approach as compared 
with that by Mello and collaborators is that it operates with an energy- dependent S*— matrix S{E). As such, it allows 
to study not only spectral correlations of different physical quantities but, in principle, also contains information 
about such features of the S"— matrix as resonances in the complex energy plane £ = E + lY. 

The notion of resonances, representing long-lived intermediate states of open system to which bound states of its 
closed counterpart are converted due to coupling to continua, is one of the most fundamental concepts in the domain 
of quantum scattering p9| . On a formal level resonances show up as poles of the scattering matrix occuring at complex 
energies £k = E]^ — ^Fj,, where E^. and are called position and widths of the resonance, correspondingly. 

The general problem of determining the domain of concentration and the distribution of poles of the S"— matrix in 
the complex plane is of fundamental interest in the general theory of scattering Jsot . Powerful numerical methods 
are available (e.g. the method of complex scaling [pT|) allowing one to extract resonance parameters for models in 
atomic and molecular physics. Whereas the issue of energy level statistics in closed chaotic systems was addressed in 
an enormous amount of papers (see [ p9|j2^ and references therein)jStatistical characteristics of resonances are much 
less studied and attracted significant attention only recently P Jl^ , [l7| , |5^ -|66| . 

In the case of weak coupling to continua individual resonances do not overlap: (P) <C A. Under these conditions 
one can use a simple first order perturbation theory to calculate resonance widths in terms of eigenfunctions of the 
closed system (see example in the Sec. HI of the present paper). Quite generally, one finds in such a procedure that 
the widths for the chaotic system with M open channels are distributed according to the so-called distribution: 

where stands for the resonance widths normalized to its mean value: — r/(P), the parameterz^ = M [v — 2M) 
for systems with preserved (broken) time reversal invariance, and V{z) in Eq.(^) stands for the Euler Gamma-function. 
The case j/ = 1 is known as the Porter-Thomas distribution pTf . It was shown to be in agreement with experimental 
data in neutron-nuclei resonances |^^, the fluorescence excitation spectrum of the NO2 molecule resonance 
dissociation of HO2 molecule [ |62[ , the diamagnetic Rydberg spectrum in lithium atom p^ ] and in microwave cavities 
1^^ . Indirectly that distribution manifests itself in fluctuations of tunneling conductance through ballistic quantum 
dots 1610. 

When the coupling to continua increases resonances start to overlap and the simple perturbative result Eq. (|^) loses 
its validity. Finally, when the coupling to continua exceeds some critical value, the so-called " trapping phenomenon" 
occurs: M very unstable states (broad resonances) are formed, whereas the rest N ~ M resonances go back to the 



real axis, i.e. become more and more narrow with increasing coupling, see [£3 59 60|] and the end of Sec. HI for a 
more detailed discussion. Such a "reorganization" of the spectrum is the most pronounced when the number of 
channels M is of the same order as the (large) number of resonances N . This range of parameters M cx: TV ^ 1 
always corresponds to the condition (P) ^ A which is just the opposite limiting case as compared with the domain 
of validity of the distribution. Under this condition one can calculate the density of resonance poles analytically 
pa,p8|. However, frequently one encounters the case of few open channels and moderately overlapping resonances 



In this situation, which is in some sense generic, one can neither rely upon the distribution 
Eq. (^ , nor use the results of (5^,^^ . The general distribution of resonance widths describing a crossover from isolated 
to overlapping resonances was found recently by the present authors for the particular case of an open chaotic system 



7l| ] . One of the main goals 
Dution, also for the case of 



with broken time-reversal invariance coupled to continua via M N equivalent channels | 
of the present paper is to give a detailed derivation and subsequent analysis of that distri 
non-equivalent channels. 

Apart from the S*— matrix elements and S*— matrix poles, the set of scattering phase shifts 9a (deflned via the S 
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matrix eigenvalues expz^a; a = 1, 2, M ) are intensively used to characterize the chaotic scattering, see |l|, [45| , [72| 



Quite recently, their statistical characteristics were studied numerically in some detail for chaotic |73 7j] as well as 
for disordered |7^,0 systems. The derivatives of phase shifts over the energy Ta = dOa/dE (we propose to call them 
" partial delay times" ) are particularly interesting and related to the mean time spent by a quantum particle in the 
interaction domain. 

The issue of the time scales associated with different stages of the quantum scattering process (e.g. tunneling, 
reflection and transmission) is quite a controversial subject which is under an intensive discussion for a long time, 
see ifrTj-^ and references therein. In particular, ambiguities arise because there is no a self-adjoint time operator in 
Hilbert space, analogous to the position operator; instead, the wave function depends on time as a parameter. 

Relegating all the essential details and derivations to section II, we just mention here that if ■0(x, t) denotes a 
wave packet at time t for a quantum particle moving in a potential C/(x) (as such satisfying the Schrodinger equation 



dt 



2m 



A + C/(x) 



'0) then the real number: 



dt / I V(x,0 P d?-> 

"'|x|<r 



may be interpreted as the total time spent by this state during its evolution inside the ball of the radius r centered 
at the origin ( we assume i — 0) normalized to unity). If il)f{-K.,t) denotes a freely evolving wave packet (i.e. 
solution of the same Schrodinger equation with C/(x) = and condition: ipf{t = — oo) ^{t = — oo)), the difference 
r(r) ~ ir[^] — tr[i^f] corresponds to the time delay inside the same ball due to scattering by the potential U{yi). The 
global time delay is defined as Td = T(r oo) and under quite general conditions (see e.g. [79| ) can be shown to be 
equal to the time-independent expectation value 

Td = j dyi-L y"dx2?A*(xi,<)rrf(xi,X2)V'(x2,0 (7) 

where Tci(xi,X2) are matrix elements of a Hermitian time delay operator Td in the position representation. This 
operator turns out to be commuting with the Hamiltonian and intimately related to the so called Wigner-Smith time 
delay matrix ||8l[| defined in terms of the S*— matrix as : 

dS^ - 

f^{E)=th—S (8) 



In particular, following the papers 82| we show in section II that the eigenvalues of the operator Td just coincide 
with the eigenvalues of fw{E). The quantum mechanical expectation value of the time delay averaged over different 
channels turns out to be equal to p^-pSl: 



^ A'/ 

a— 1 

where the bar denotes the averaging over the energy spectrum of the packet. This shows the relation between the 
phase shift derivatives and mean time delay mentioned above. 

A quite detailed analysis of the time-delay problem was given in the context of nuclear physics by Lyuboshits [^3|j84[ 
and other authors |8^. In particular, for a wave packet of arbitrary form Lyuboshits [Q suggested a concept of the 
probability distribution of its time delay. His definition is based on the interpretation of the quantity P{t) — J^^y \ 
'(/'(x, t) p d'^x as the quantum mechanical probability to be found within the volume V at instant t. Then the time 
derivative dP/dt can be used to define the distribution of times spent inside the volume V. A general and illuminating 
discussion of the time evolution properties of wave packets in a generic chaotic systems can be found in [ ^2| . 

On the other hand, the existence of the Hermitian time-delay operator Td in Hilbert space suggests an alternative 
definition of the time delay statistics by the natural requirement that tP = {'^{t) \ T^ \ ^(t)) for any wave packet 
^(t). Then the problem is reduced basically to study the statistical properties of the Wigner-Smith time delay matrix 
tUE). 

The chaotic scattering makes the Wigner-Smith time delay matrix (in particular, the quantity j^Trr^ which is 
called just Wigner time delay) to be a strongly fiuctuating function of the energy E as well as of any external parameter 
X. From this point of view we can speak about distributions and correlation functions of these quantities. Similarly, 
the distribution of partial delay times can be used to characterize variations of time scales associated with the chaotic 
scattering process. Various statistical aspects of time evolution of the chaotic quantum systems were studied earlier 
in some details in [|2|J8|J8|]. 
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Being an important characteristic of the scattering process, the statistics of phase shifts and their derivatives deserve 
to be investigated in more detail. Additional interest to the problem attaches the fact of relevance of Wigner time 
delay in condensed matter physics. Indeed, in a series of papers by Biittiker and collaborators the Wigner time delay 
was related to the frequency-dependent response of mesoscopic capacitors l8^-|9^ . A more detailed discussion of this 
issue can be found in the last section of the present paper. 

More general parametric derivatives of the scattering phase shifts can also be related to some observable quantities. 
As a particular example we mention the relation between the persistent currents and the derivative of the total phase 
shift over the magnetic flux derived by Akkermans et al. These authors considered "open mesoscopic networks": 
two-dimensional systems of conducting loops coupled to infinitely long ideal leads (waveguides). The loops can encircle 
a flux tube with flux <j). The expectation value of the persistent current around the flux tube in the state | Vf) is 
/(4',(/)) = — (vlf I dH/dip I ^). For the case of a closed system (i.e. when loops are disconnected from the leads) 
each discrete level i?„((/)) carries the current —dEn/dcf). When the system is open , it turns out that the differential 
contribution of the scattering states at energy E to the persistent current can be expressed in terms of phase shift 
derivatives as [0: 



dI{E, 



J_d_ 

2m dd) 



lndetS'(^;,. 



dE 



Y^e{E,. 



dE 



(10) 



In the present paper we give quite a detailed analysis of the statistical properties of scattering phase shifts and 
their derivatives for generic chaotic scattering in a system with broken time-reversal invariance. The extension of our 
results to other symmetry classes as well as to the crossover regimes will be published elsewhere |^] . We find it to 
be most informative to concentrate our attention on the so-called /C— matrix related to the scattering matrix as: 



5* = 



I + iK 



(11) 



This equation shows that eigenphases Oa{E, X), a = 1, Af considered modulo 27r are determined in a unique way 
by the eigenvalues Za{E,X) of the iiT— matrix , where we have indicated explicitly both the energy dependence and 
dependence on an external parameter X. First of all, we calculate explicitly the correlation function of densities of the 
eigenvalues of the X— matrix at two different energies E ±0. and parameter values X ± 5X . When 17 = 5X = this 
correlation function turns out to be the same as the pair correlation function following from the Poisson's kernel. This 
fact confirms the expected equivalence of the minimum information ap pro ach and Hamiltonian approach for the case 
of fixed energy, and as such extends the earlier studies on that subject |l6| . From that moment wc concentrate on the 
statistics of delay times and parametric derivatives of phase shifts. First, we derive and analyze the general expression 
for the distribution of "partial delay times" = dOa/dE (here and henceforth we frequently put ?i = 1) as well as 
derivatives dOa/dX . This distribution, being an interesting characteristic of the chaotic scattering by itself, also allows 
us to detect the qualitative features of the Wigner time delay distribution. In particular, for the one-channel system 
the partial delay time is exactly the same as the Wigner time delay. After that, we derive the parametric correlation 
function of the Wigner time delay and show some interesting correspondences with the results of the semiclassical 
approach. A short account of our results was published earlier | |7l| , p3[ |. 

The organization of the paper is as follows. Section II is meant to be a kind of introduction to the random matrix 
method of the descr iption of an open chaotic quantum system. It is based mainly on the original papers by other 
authors 
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94| , |9q |. Considering a particular generic example we discuss the main ingredients of the model and 
present a quite detailed discussion of the time-delay operator and other quantities characterizing time evolution in 
such systems. 

In Section III we derive and analyze the density of the S-matrix poles in the complex plane. Section IV is devoted 
to the statistics of eigenvalues of the X— matrix, phase shifts and their derivatives and analyses different statistical 
aspects of the Wigner time delay. It contains also a kind of semiclassical analysis of the parametric correlations of 
Wigner time delays. Concluding remarks and a discussion of the potential experimental relevance of the obtained 
results can be found in the final Section V. 



II. SCATTERING PROBLEM FOR RANDOM-MATRIX HAMILTONIANS 



A. General description of the model. 



A model which is most appropriate for incorporating random matrix ideas for describing the phenomenon of 
quantum chaotic scattering was discussed in great details in the works by Verbaarschot, Weidenmiiller and Zirnbauer 
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| p2[ and Lewenkopf and Weidenmiiller . A general construction actually goes back to works by Feshbach and 
is based on the theory of quantum scattering as formulated in the book by Mahaux and Weidenmiiller on the 
theory of nuclear reactions, or in the book by Levine about quantum scattering in molecules. Let us also mention 
a very general and profound analysis of open quantum and classical wave systems performed in the book by Livsic 
psf . In order to make the present paper self-contained we give below a list of the main ingredients of the model. 
Since the general construction is a rather abstract one, we find it instructive to illustrate it by presenting the explicit 
derivation of the expression for the scattering matrix and its subsequent analysis for a generic system of much interest: 
scattering of a quantum particle moving in a perfect lead of width d in contact with a chaotic region, simulated by 
a random matrix Hamiltonian. Such a derivation follows that of the papers |^,^ (see also the book |Q and the 
works by Pavlov and collaborators ]40|]). 

When dealing with a generic scattering problem it is natural to assume that the scattering event is always confined 
inside a compact part of the available space which is called the " interaction region" . Outside this region interaction is 
absent and fragments exhibit a free motion characterized (apart from the total energy E) by a set of quantum numbers 
describing the internal quantum state of each fragment. As such, these quantum numbers specify the states in which 
particles can be found long before and long after the scattering takes place - the so-called " channels of reaction" . For 
example, in nuclear and molecular physics different channels are marked by relative angular momentum and spins 
of colliding particles. In the particular example considered in much detail below the motion of a particle along an 
infinite lead of width d is quantized along the transverse direction, the different transverse modes being different 
channels. Assuming that exactly M channels at given energy E are "open" (i.e. allow for an unbounded motion of 
the particles), we associate with the channel region a continuous set of functions \ a, E); a = 1, M normalized as 
(a. El I b, E2) = 5ab5{Ei — £2). An analogous, but discrete set of orthogonal states \ n); n = 1,2, N is associated 
with the compact interaction region. 

In the absence of interaction between the states in channels and the internal states the Hamiltonian of the system 
has obviously the form: 

Tio = E I iH,n)un {m\+Y, [dE\a, E)E{a, E \ (12) 

where the integration goes over the energy region where the given channel a is open. The model is simplified by 
neglecting any direct coupling between different channels; hence the corresponding term in the Eq.([l^) is diagonal 
in a. The first term describes the Hamiltonian of the "closed" chaotic system possessing iV ^ 1 bound states, 
which are eigenstates of In the spirit of the Random Matrix universality conjecture, we simulate this part of 

the Hamiltonian by taking Hin to be a Gaussian random N x N matrix. The number N is considered to be large: 
N ^ 1. To describe the interaction between channels and bound states (converting bound states into resonances) one 
adds to the Hamiltonian Eq. (O) the " interaction term" : 



T> = ^ M / dEWia{a,E \ -|-herm.conj 

7 ^ \ 



(13) 



In any practical implementation of such a procedure one should make sure that the total Hamiltonian Ti. = Ti.Q + V 
is self-adjoint. This point is not at all trivial (see the example below). The general way of self-adjoint matching of 
Hamiltonians with internal structure and those describing motion in external scattering channels was suggested by 
Pavlov and developed by Pavlov and collaborators pQ] . A particularly convenient formulation of the Pavlov's method 



suggested by Makarov |41| was applied to the problem of chaotic scattering in the recent papers 

After the self-adjoint matching is done one can employ standard methods in scattering theory "(see the mentioned 
books [t36y37| ]) in order to write down the Lipmann-Schwinger equation for the in-and outgoing scattered waves and 



find an explicit expression for the scattering matrix. 



B. From Random Matrix Hamiltonian to Scattering Matrix. 

Instead of demonstrating such a formal derivation within a quite abstract "projection formalism" (see e.g. [p5| , ^ ) 
we find it to be more illuminating to show how to derive the 5— matrix in an alternative way |^,^. To do this 
let us confine ourself to a particular generic example of a scattering system depicted schematically in Fig.l.: a two 
dimensional cavity of irregular shape with impenetrable walls coupled to an infinite waveguide (lead) of width d. 
Let us mention, that it is one of the favorite models for the study of generic features of chaotic scattering, both 
theoretically (see e.g. the "frying pan" model in |Q) and experimentally |^,|l^. The propagation of a quantum 
particle inside the lead is described by the Schrodinger equation: 
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+ = *(x,, = ±./2).0 (14) 

whose general solution can be represented as a M— component vector: ^ = (^'i(x, y), .... 'i'mix, y))^ with components 
^a{x,y) = i)a{x)(j)a{y) , whcrc 

'^-(^) = ,0 .2;\ u/2 [^ae-^'^"^ + Bae^'="1 ; UV) = (2/rf)^/^ sin [(^) (y + d/2)l (15) 
for X > 0; I 2/ |< (i/2; a = 1, 2, M, with the number M of open channels at the energy E — being equal to the 



largest integer less or equal to and the wave vector ka being equal to ka = fc^ — (^)^ , so that —-^ipa = k'^4'a- 

The running waves are properly normalized to energy J— functions; unitarity of the 5— matrix to be introduced later 
is related to the conservation of the probability flux. 

The situation of the waveguide disconnected from the cavity we describe by the boundary conditions: dipa/dx \x=o= 
0. This means that the particle in each channel is just reflected back: Aa = Ba. The corresponding 5— matrix relating 
the vectors of incoming A — {Ai, ...,Am)^ and outgoing B = {Bi, ...,Bi\j)^ amplitudes: B = SA is just the unity 
M X M matrix: S — 1. The role of the vector | a, E) of general construction (see Eqs.( p^p^ )) is played by the vector 
^ corresponding to the particular choice of amplitudes of incoming waves: Aa = 1; Ah^a — 0. The Hamiltonian of 
the particle motion inside the cavity is simulated by the N x N random Hermitian matrix ■ Correspondingly, the 
"internal" wave function is represented by the A^— component vector u — {ui, ...,ui\[)'^. The wave functions of the 

scattering system as a whole (" cavity attached to the lead" ) are therefore vectors ^ = ^ ^ ^ from the Hilbert space 

X^(i?+, C^"^) ® supphed with the scalar product: 



1/2 



/d/2 noo 
dy / da;*|*2 
-d/2 Jo 



'-d/2 

Let us define the Hamiltonian operator H of the system as a whole acting in that Hilbert space as: 



, u \ / H„^u + J^^Ldy dxyV{x,y)^{x,y) \ , , 



where Ti-ch is the operator diagonal in the channel space: 



and yV(x, y) and V{x, y) are N x M and M x N rectangular matrices describing a coup ling be tween two parts of the 
Hilbert space. Let us assume for simpHcity that the coupHng is local along the waveguide ||94| , |95|] : >V(a;, y) — 5{x)W{y), 

so that f^^f/^dy dxyV{x,y)^{x,y) = J'^^f/^dyW{y)'if \x=o^ w^|J{x = 0), where Wia = J'l^f/2dyy^ia{y)<Pa{y) i = 
l,...,iV, a — 1,...,M and il){x) = (il;i{x), ...,iPm{x))'^ . On the other hand, we have to put V{x,y) = in order to be 
consistent with the locality of the coupling and to stay in the space spanned by the vectors 
The operator Ti defined in such a way is not, in general, a self-adjoint one. Indeed, 

(H*l, *2) = Ui^i?Lu2 + (Tich*!, *2) +4'lix = 0WU2 (17) 

and 

(*1,H*2) = UitiJ,„U2 + (*i,7T;,ft*2) + ulu>V2(.T = 0) 



From the definition of the operator Hch and that of the scalar product (vPi, '$'2) one can easily find after using partial 
integration and the fact Hj^ — Hin that: 

(H*l,*2) - (*1,W*2) = (18) 

- (^V'2)} U=o = 0)7«tu2 - n\wMx - 0) 
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In order to have a self-adjoint Hamiltonian operator H. one has to impose some appropriate boundary conditions 
at the point cc = ensuring that the expression above is vanishing [p6| . The most obvious (however, not the most 
general) choice is: 



d 



(19) 



On the other hand, the solution of the Schrodinger equation for the whole system: H.^ = E — (we call 



2m 

these solutions the " scattering states" and denote them $ e henceforth) leads immediately to the relation: 

-1 



E - H,n) wi^ix = 0) 
which together with Eq.([T9|) yields the following equation for the vector '4'{x): 



h id 



(20) 



(21) 



This equation allows us to find easily the unitary scattering matrix: 

-1 



I-iK 



X 



I + iK 



; K = wW^ 



E — Hi. 



-W 



(22) 



where W — ^ 2m / t:Ti^ wd\a.g{k^ /cj,/^^), / is the unity matrix of the corresponding dimension. Often we treat 

diagonal matrices proportional to / simply as numbers. Sometimes, we indicate the dimension as index, here it is 

The expression Eq.(|22|) can be also rewritten in another form, frequently used in applications. To this end we write: 



I-iK 



I + iK 



I + iK - 2iK 



I + iK 



= 1 — 2« ( / + iK ) K and use the identity 



i + inW^- 



E — Hi, 



-W 



= ^(-«7 

E - H,n f^^f^ 



1 



I + iix- 



1 



1 



E — Hi: 



-W = 



E — Hij 
1 



-W 



k+l 



(23) 



E — Hin 

which means that the scattering matrix can be written in the form: 

S = i-2i7rW''- ^ 



E — Hi, 



tWW^ 



-W 



E~H. 



-W 



(24) 



where the non Hermitian effective Hamiltonian Tie f is given by Hef = H„i - iT and f = nWWK 

The expression for the scattering matrix of the form Eq. (|24) a ppears generally when one describes an open quan- 
tum system decaying into several open channels, see e.g. ]38|, 39| . In particular, it can be derived from the general 
Hamiltonian Eq. ( p^ , p^ under the assumption that the elements of the matrix W are energy-independent. This was 
just a starting point in the approach by Weidenmiiller and collaborators E^,|2^. In the derivation above the matrix 



We, however, will be mostly interested 



W does depend on the energy E via the parameters ka — 

in the fluctuation properties of the energy-dependent S'-matrix. The typical energy scale of such fluctuations is given 
by the mean level spacing A - a typical separation between the adjacent eigenvalues of the matrix Hin. Far from 
the thresholds, as long as A is negligible in comparison with the difference between the adjacent threshold energies 
1^(^M ~ ^M-i) ~ '^^m^ (if) ^® '"^'^ safely neglect the energy dependence of the matrix W. In view of A c>c 1/A^ 



the latter requirement is always satisfied in the limit iV 3> 1 which is the only case studied in the present paper 
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C. Time Evolution of the Wave Packets: Staying Probability and Time Delay Operator for Open Quantum 

Chaotic Systems 



Before describing the time evolution of wave packets let us note that any particular scattering state = 
{u{E) , ^ (E))'^ is uniquely specified by the set of incoming amplitudes A = {Ai{E), Am{E))'^ . Being the eigen- 
functions of the Hermitian Hamiltonian Ti the scattering states must be orthogonal. Below we verify by a direct 
calculation the validity of the orthogonality condition: 



- U^(^i)u(^2) + = S{Ei - S2)AIA2 

Such a calculation allows one to derive some helpful relations that are used later on. 
By using the identities 



(25) 



1 

2^ 



dxe 



1 



6iu)- 



1 



1 



2mu ^ka{EME2) 



S {kjE,) - ka{E2)) ^ SiE, - E2) 



and exploiting the definition of the S'-matrix: B(i?) — S{E)A{E) and its unitarity one easily finds: 

1 



i^E^,'^E,)^S{Ei-E2 



) (aI A2) 



27ri{Ei - E2) 



A\ (^S\Ei)S{E2 



DA. 



(26) 



Here we assumed that both energies Ei, E2 are far from thresholds and close to one another, so that effectively we can 
put ka{Ei) = ka{E2) in the expression above everywhere (this is consistent with neglecting the energy dependence of 
W as discussed earlier), except in the denominator, where {Ei — E2) oc (fc^(i?i) — fc„(£'2)). 
Let us now use the relation Eq.(C^ 



rewritten in the following form: 
1 1 



2 ^ - i/„ 



W [I + S{E)] A 



,t 



Hence: 

Now we use the identity: 



-a\ [i + s\ES) 



1 



1 



1 



E^ — Hin Eo 



Hi, 



-W = 



El — Hin E2 — Hi-, 



-w 



{l + SiE2)) 



El 



W = 



El — Hin E2 — Hi, 

—J——(k{Ei)~K{E2 

TT{E2 - El) \ 



(27) 



(28) 



(29) 



where we exploited the definition of the K matrix, see Eq.(|22|) neglecting the energy dependence of the matrix W . 
Relation Eq.(|2^) between S'-matrix and i^T— matrix can be written also as K{E) = — S{E)j ^/ + S{E)^^ 

(E) . Substituting this relation into Eq. (|2^) we use it to reduce Eq. (^8|) to the following final form 

1 



2TTi{El-E2)' 



A\(s\Ei)S{E2)-i) A 



(30) 



We see that when combined together the Eqs.(|26|) and ( pOj) produce exactly the orthogonality condition Eq.(^5|). 

In particular, this orthogonality condition allows us to use the scattering states ^^-^ corresponding to the choice of 

incoming amplitudes Aa — 1, A^^a = as a convenient basis in the full Hilbert space. Denoting =| #^'') one 
can write down the total Hamitonian Ti. as : 



» M 

J dEY,\^^^'>)E{^ 



(a) 
E 



(31) 



Now we are prepared to answer the following question: given a wave packet | $(t)) = (u(i), ^{t))'^ which evolves 
according the Schrodinger equation ih{d/dt) \ ^{t)) = TL \ ^{t)), how to express in terms of the scattering matrix 
the probability for the corresponding particle to be found inside the "chaotic" domain a; < at instant t. 
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According to the rules of quantum mechanics this probability is just given by P{t) = u^{t)u{t). Let us expand the 
wave packet over the scattering states | = X^fli \ ^b'') 

I = J dEf{E) I *^)exp-^£;; J dE \ f{E) 1 (32) 

where the coefficients f{E) determine the initial form of the wave packet and we assume: A^^A = 1 so that (<i>(t) | 
<i>(t)) — 1. This immediately gives us the desired expression: 

P{t) = J dE,dE2f*iEi)f{E2WEUE, exp-|(Si - E2) = (33) 

ab 

where we made use of the Eq.(|30|). The mean time spent in the interaction region (i.e. the mean time delay r^) can 
be found integrating this expression over the time. This operation produces a S— functional factor S{Ei — E2) in the 
integrand which finally gives: 

/C30 
p{t)dt = ^:A^a^ (34) 

-°° a,b 

where 

r{E) = ^n^^SiE) = -^hS^E)^ (35) 

is the Wigner-Smith time delay matrix [p7| and the bar stands for the energy averaging determined by the wave 
packet spectrum: {■■■) — J dE{...) \ f{E) p. If the particle comes only via a particular channel a the scattering states 

I ^e) coincide with the basis states | ^^■*) and the corresponding time rjf ^ (which is natural to call the "delay time 
for the channel a") is given by the spectral average of the diagonal element faa{E). Then the delay time averaged 
over all channels is given by the spectral average of the Wigner-Smith time delay: ■^Trf(£^). 

The time derivative —dP/dt is a current out of the chaotic region. Assuming that the spectral function f{E) varies 
with i? on a much larger scale than the mean level spacing A ( the latter scale is typical for variations of 5*— matrix 
elements) we can put f{Ei) w /(-E2) in the expression Eq. (|33|) This results in the following expression: 



d_ 



P{t) ^ S{t) ^ Pit); p{t)^^Jdee^^'/'^Y.^aAb{{sHE + e)S{E))J (36) 

ab 



This expression can be interpreted as follows [ p4| , p2[ . In our approximation the part of the Hilbert space corre- 
sponding to the chaotic region is not populated at t < 0. At t — the wave packet reaches the chaotic region and 
populates its states instantly. This fact is described by the d{t) term in the expression Eq.(^). Then the function 
p{t) has a meaning of the distribution of duration of stay inside the chaotic region. This was just the reason to call 
p{t) the distribution of time delays p^ . On the other hand, according to conventional rules of quantum mechanics in 
order to speak about the probability distribution of some observable one should be able to find a Hermitian operator 
in Hilbert space generating all the moments of that observable as expectation values of the integer powers of this 
operator. Some important insights in the issue of constructing such an operator can be found in |79| . To this end let 
us consider the following time delay operator constructed in terms of Wigner-Smith time-delay matrix Eq.(Bq) as: 



dEY,\^^E^)ir)abiE){^'^^^\dE (37) 

2fc 



The Hermiticity of this operator follows from that of the Wigner-Smith matrix. It commutes with the Hamiltonian, 
Eq. (^l|) , but is not at the same time diagonal due to the degeneracy of Ti. It is evident that for any wave packet 
I ^{tf) the mean time delay given in Eq.(p4) is just the time independent expectation value: 

(t) ^ r, - (*(t) I fw I m) (38) 
Then it is natural to define the higher moments of the time delay as : 
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a.b 



(39) 



This should be contrasted with the moments of the distribution function p{t): 



I dtt^pd) = {^/nY ^ AlA, ( ^^^^S{E) j (40) 

a.,b \ J ah 

We see that only the first moment of this distribution coincides with that given by Eq.(^9|), all other being dif- 
ferent. This particular example shows certain ambiguity in definition of delay time statistics. In the present paper 
we concentrate on statistics of Wigner-Smith time delays and related quantities: energy derivatives of S*— matrix 
eigenphases. 



D. S-matrix characteristics: poles, eigenphases and delay times 



The expression (24) forms the basis for extracting the statistics of scattering poles (resonances), which are merely 
the complex eigenvalues of the non-Hermitian effective Hamiltonian Tig/ = Hin — inWW^ . At the same time, the 
expression (|2|) turns out to be a more convenient starting point for studying statistics of scattering phase shifts and 
delay times. Indeed, it is evident that scattering phase shifts 9a', a — 1, 2, M [defined via the S-matrix eigenvalues 
exp i9a] are determined by the eigenvalues Za{E, X) of the matrix ttK in view of the relation: 9a — —2 arctan Za- Here 
we indicated explicitly the dependence of the eigenvalues Za on the energy E and an external parameter X originating 
from the corresponding dependence of the Hamiltonian on the parameter: Hm = Hin(X). It is therefore convenient 
to characterize the statistics of phase shifts via the spectral density: 



M 

PE,x{z)^—Y.^{z-Za{E,X)) (41) 

a=l 

Actually, the relation 9a = —2 arctan determines the phase shifts modulo 27r only. It is easy to understand that 
every time the energy E coincides with one of the eigenvalues En{X), n — 1, 2, TV of the Hermitian Hamiltonian 
Hin{X) one (and only one) of the scattering phase shifts crosses the value 27r x integer. Indeed, assuming that the 
eigenvalues En are generically not degenerate one can write the matrix elements of the matrix K in the vicinity of 
E = En as K{E — > En) = W*aWnbi where Wnb are matrix elements of the coupling matrix W in the basis of 

eigenstates | n) of the Hamiltonian matrix Hin and (l^^)^„ ~ ^na- We see immediately that K{E ^ En) has only 
one eigenvalue divergent at _E i?„ which is given by z{E) = -gz^ X^a I the corresponding eigenvector 

being v = (W^i, W^2: ••■j l^nA/)"^- The phase shift corresponding to the infinite value of z{E) must be an integer of 
27r. At the same time all other eigenvalues of K{E — > En) are exactly zero in that approximation with corresponding 
eigenvectors belonging to the M—1 dimensional space orthogonal to v. This fact just means that M—1 corresponding 
eigenvalues of the exact matrix K{E) stay finite in the vicinity of En- Introducing the exact density of states for the 
closed chaotic system: vx{E) = jj J2n ^ ~ En{X)) and fixing the phase shift value a,t E — — oo to be zero, we 
conclude that: 

M „E 

^9a^ 2ttN / dwx (m) - 2 ^ arctan Za (42) 

Here arctan Zq means the principal branch: | arctan z^ |< 7r/2. As function of energy E this expression is continuous 
and monotonically increasing with E. The first term is proportional to the level staircase and we can forget it modulo 
2tt. 

We will use the relation ( p2| ) later on in order to determine the correlations of the Wigner delay times Tw{E) = 
{d/dE)jjJ2a^a. = —{i/M){d/dE)ln'DetS{E), which , of course, are positive. For the latter quantity we also can 
find an independent representation by noticing that: 

In DetS{E) = In '^J - In 1. , . (43) 

I)et{I + iK) Bet{I + iTT^-^VVVVf) 
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where we have made use of the identity 



Det(/ - UV) = Det(/ - VU) (44) 
vahd for arbitrary (also rectangular) matrices U,V. From Eq.(|43|) we immediately obtain the simple relation: 

rUE) ^ -(2/Af)Im Tr {e - iJ,„. + ^.WW^) ^ A f Jf (45) 

which, in particular, shows an intimate relation between the statistics of Wigner time delay and that of 5-matrix 
poles. 

To study 5— matrix characteristics within the framework of the stochastic approach, one should specify the proper- 
ties of the amplitudes Wia, which couple the internal chaotic motion to M open channels. For the sake of simplicity 
one can restrict the consideration to the case when the S— matrix is diagonal after averaging: {Sab) = Sab{Saa)- Such 
a choice is related with the mentioned absence of direct coupling between the channels |^,|2^ and can be ensured if 
one considers fixed amplitudes Wia satisfying the so-called orthogonality relations Q : 



^ W:aW,b ^ -JaSab (46) 
TT 



An alternative way to ensure the diagonality of the average 5— matrix is to consider the amplitudes Wia to be 
independent Gaussian random variables p3| , ^ , ^7|j5^ ] : 

7a 



{W,a) = 0; {w:aW,b) = ^SabS,, (47) 



One can show, following the papers that both choices lead to the same results as long as the number M of 

open channels is negligible in comparison with the number of bound states TV. Since this case is the only considered 
in the present paper, we restrict ourselves to the condition Eq.(^) henceforth. Provided the orthogonality condition 
( ^ ) is fulfilled, one can show that the diagonal elements of the S*— matrix are given by the following expression: 



where g{E) = iE/2 -|- (1 — E'^/4)-'^/^. We do not give here the derivation of Eq.(48) referring the interested reader to 
the paper for more details. However, we mention that for the one-channel case M — 1 relation Eq.(^) follows 
directly from the distribution of the phase shift 9 to be found in Sec. IV. The strength of coupling to continua is 
convenient to be characterized via the transmission coefficients Ta = 1 — |(S'aa)P that are given for the present case 
by the following expression: 



2 



(49) 



The quantity Ta measures the part of the flux in channel a that spends substantial part of the time in the interaction 
region |^,^ . This interpretation follows from the fact that the energy averaged S— matrix (equal to the ensemble 
average (S) by ergodicity requirement) describes a short time scattering ("direct response", see the Introduction). 
Let us also note that frequently we find it to be more convenient to use the parameters 

5a - |- - 1 (50) 

rather than the "transmission coefficients" Ta . 

Naively, one could suspect that the larger is the parameter ja, the larger is the part of the flux effectively penetrating 
the chaotic region. However, we see that this is not the case: both limits 7a ^ and 7a ^ oo equally correspond 
to the weak effective coupling regime <C 1 whereas the strongest coupling (at fixed energy E ) corresponds to the 
value 7a = 1. The maximal possible coupling corresponding to the upper bound Ta = 1 is achieved in the present 
model for an energy interval in the vicinity of the center E ^ 0. 

This feature is a purely quantum effect and is not surprising any longer, if one remembers the simplest textbook 
example of a quantum particle scattered on a one-dimensional " potential step" : V{x) = for 2: < and V{x) = V 
for X > 0. The transmission coefficient for such a problem is given by: T = 1— | 5 p= 4kK/{K + fc)^, where 
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k = {2mEf/^/n; K = [2m{E - V)^/^ /Ti and S > stands for the energy of incoming particles. Similarly to 
the case above the transmission is very small both for the system "almost closed" classically: E — V <^ E when 
7o = {K/ky/'^ ^ 1 as well as for systems "very open classically" V < 0,\ V \^ E when 70 3> 1, the "perfect" 
transmission T — I being possible for the only case = when 70 = 1. This simple example is of course just to 
remind us of the effect known in radiophysics as " impedance mismatch" : the wave is always reflected back at the 
point of contact of two different waveguides, unless special boundary conditions are ensured. 



III. SCATTERING POLES IN COMPLEX PLANE: DISTRIBUTION OF RESONANCE WIDTHS. 
A. Resonances as eigenvalues of a non-Hermitian Hamiltonian. 

We are interested in determining the average two-dimensional density 

of complex eigenvalues £j = Ej + iYj, j = 1,2, ...,N of a non-Hermitian effective Hamiltonian Tie/ = Hm — zF, 
see Eq.(p^. According to the general discussion presented in the Introduction we use a, N x N random matrix Hin 
from the Gaussian Unitary Ensemble to model the Hamiltonian of a closed chaotic system with broken time-reversal 
symmetry. The entries Tij of the matrix F are expressed in terms of the channel amplitudes Wia, a = 1,2, ...,M 
as Tij — TT^^ WiaW*^. Before presenting a general theory it is instructive to consider the important limiting case of 
an extremely weak coupling when we expect that the resonances are so narrow that their widths are much smaller 
in comparison with the mean separation between the unperturbed levels. Under these conditions a simple first order 
perturbation theory is adequate. Using the notation | n) for the eigenvector of Hi„ corresponding to the real eigenvalue 
En of the closed system: Hi„ \ n) — £"„ | n) , one can estimate the shift of the eigenvalues into the complex plane as 

N 

- y„ = (n I f I n) = ^ (a/ I f I ak){n \ ai){ak \ n) (52) 

where | a^); k — l,2,...,iV is an arbitrary chosen basis of orthonormal vectors. The matrix F can be easily 
diagonalized and shown to have exactly M non-zero eigenvalues 7a = "'X^i ^L^ia- Choosing | a^) to be eigenbasis 
of the matrix f we therefore have: Yn = — X^aLi ^a,{n \ oia){aa \ n). Now we use the well-known fact that different 
components (n | aa) = Ua+iVa of eigenvectors of the GUE matrices in an arbitrary basis can be treated as independent 
complex variables, their real and imaginary partsbeing independently distributed according to the Gaussian law with 
the variances = = 2^ . This fact allows one to calculate the distribution of Yn easily. Considering for simplicity 
the case of all equivalent channels: 7a = 7 for any channel, we get: 



V{Y) = d{Y - Yn) = 



^ikY dk 



2tt 



N 

dudv — exp ^{N — ijk){u^ + u^) 

TT 



(53) 



dk 1 _ N^' \N\Y\ 

2tt'' {1-ikj/Nyi 7^^F(M) ' ' '^''P" 
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for F < and zero otherwise. 

We arrive at the well-known result: the widths of resonances for a slightly open chaotic system is given by the 
so-called distribution |^l[]. Actually, the same form is applicable also for M-channel open systems with preserved 
time-reversal symmetry, provided one changes M — > M/2. The latter distribution for M = 1 is known as the 
Porter-Thomas distibution. 

When the coupling to continua increases some resonances start to overlap and the simple perturbation theory loses 
its validity. 

A general method for calculating the eigenvalue density for non-Hermitian random matrices was proposed by 
Sommers et al. [Q. These authors suggested to recover the density p{E,Y) from the "potential function" 

- <f(i?,r) = 1 (inDet {£ - Hef}^ (£ - Hef)) = (54) 
( 1 5: In I £ - I A ^ (1 5: In [(£; - E,f + [Y - Y,f] ) (55) 
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To show that this is indeed possible it is convenient to regularize the logarithm in eq.(^5|) first and consider [5^ : 

- HE, Y),^(^Y.^n [{E - E,)' + {¥- Y,f + e^] ) (56) 



N 

3 



Then one notices that the function: 



piE,Y), ^ ~^{dl + d^)<i>iE,Y), = -r^y ^ (57) 



produces the required two-dimensional density, Eq. (|5l|) , when e — )■ 0. Indeed, for an arbitrary continuous function 
f{E, Y) one has: 

lim- / / dEdYf(E,Y) ^= (58) 

^-"^i-ooi-oo '[(E-Ej)^ + {Y-Yjy+e^f 

hm - / / dudvfieu + E,,ev + Y,)— -— ^ = fiE„Y,) 

in agreement with the S- functional property. 

In fact, the expressions Eqs.( p6| , p7| ) show that p{E, Y) can be considered as a two-dimensional density of fictitious 
pointlike "electric charges" the function ^{E,Y) playin g the role of the electrostatic potential for such a system 

and eq. (|7|) being the corresponding Poisson equation fl^^Q . 

Actually, it turns out to be more convenient to use a slightly different regularization, as it has been actually done 
in HQ 

^E,Y; k) ^-^(inDet [(£ - He/)^ (f - He/) + k'/w] ) (59) 

performing the limiting procedure k — > at the very end. For the so-called normal matrices (whose Hermitian 
conjugate Ti.^ commutes with Ti.) the regularized potential Eq.(^9[) coincides with that defined in the Eq.(|5^). It is a 
less trivial fact that one recovers the two-dimensional density of complex eigenvalues from the potential ^{E,Y;k) 
also in a general case of nonnormal matrices Ti.. We show in the Appendix ^ that it is indeed the case: the density 
obtained from Eq.([5S|) by Poisson's equation is positive and goes to a sum of (5— functions with weight 1/N near the 
eigenvalues of He/ ■ 

The main technical problem is to perform the averaging of the logarithmic potential ^{E,Y; k) over the random 
matrices from the corresponding Gaussian ensemble, eq.(^. To perform such an averaging the authors of the papers 
p5| , p8[ employed the famous, but somewhat problematic "replica trick". This procedure amounts to averaging the 
corresponding determinant raised to an arbitrary positive integer power n, the average logarithm being recovered as 
a result of the limiting procedure n — > 0. In general, however, the analytical continuation rt ^ is not unique. In 
particular, it is known that the replica trick fails to reproduce correctly the correlation function of densities of real 
eigenvalues of large Hermitian matrices at two points E ±n of the spectrum [Q. Rather, it succeeds in giving the 
correct behavior of that correlation function at the scale Q large in comparison with the typical separation between 
neighboring eigenvalues, known as the mean level spacing A. As is shown below, formally the calculation of the 
mean eigenvalue density in the complex plane is very similar to the calculation of a correlation function of eigenvalue 
densities on the real axis, with the role of Q played by the variable iY measuring the distance from the real axis. We 
immediately conclude, that the replica trick must fail when we are interested in eigenvalues situated sufficiently close 
to the real axis: F ~ A. 

To this end it is necessary to mention, that the non-Hermiticity of the matrices considered in the papers [p5 



was, in a sense, quite strong: the probability for an eigenvalue to be situated at the distance Y A (x 1/N vanished 
in the limit N ^ oo. Under those conditions it is not surprising that the replica trick succeeded in producing the 
correct result, the fact verified both by independent methods: by a variant of the supersymmetry method |^8| and by 
direct numerical computations [ p5| , p8[ . 

The situation described above is drastically different from that we expect to happen at the model under the present 
consideration. Indeed, it is known, that when the number M of open channels is small in comparison with the number 
N of relevant resonances, the majority — M of resonances are rather "narrow" and the corresponding poles are 
situated close to the real axis [ |5^ , |55| , p9| , p6[ . Under such a situation one has to discard the replica trick and to seek 
for a more reliable procedure. Fortunately, the authors of the mentioned paper |99| showed how to calculate the 
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two-point correlation function for real eigenvalues correctly by exploiting the method pion eere d by Efetov ]10C( | in the 
theory of disordered solids. A pedagogical introduction to the method can be found in |101]. In the present paper 
we adjust this procedure for finding the density of scattering poles in the complex plane for the few-channel case |7|] 
(for the many-channel case M oc N this density has already been calculated by N.Lehmann et al. [^). 

Instead of working directly with the regularized potential ^{E,Y, k), see Eq.(^9|), in terms of which the two- 
dimensional density p{E, Y) is expressed as: 



we prefer to consider the related function: 

1 



\im—id% + dl.ME,Y,K) 
K^o 47r 



N 



In Det 



E + n- Hi, 



iT] E-n- H,,r 



N 



(60) 



(61) 



It is evident, that the potential ^{E, Y, k) can be obtained from the function $(£', —iil; k) by analytical continuation 
—iil — > Y. As long as k is finite, there is a region extending from positive to negative Y where the function $ is 
analytic in y = — if2. Actually, this continuation is more convenient to perform directly on the level of densities, i.e. 
first to calculate the auxilliary function |l02[ : 



and to restore the true two-dimensional density p{E, Y) letting — — > Y first and then k — > 0: 

p{E,Y) = Yimp{E,-i^l = Y),. 

To this end let us consider the generating function: 
Z(£;,r!;£;b,r!,;K) 



Dot 




[E + n- H,n + if) 


(^E-n- H,n - ?:f ) 






Det 


[( 


'Sb + Sib - H,n + if ^ 


(Eb - fife - - 


+ K^iN 



(62) 



(63) 



(64) 



in terms of which the function p{E, — iSl)^ is expressed as follows: 



p{E, -iri)„ = — 
An 



— hm — 

dE (Ei,s2t)^{E,n) dE 



— hm — 



(Z{E,n;Eb,nb;K)), 



(65) 



The determinant in the denominator of expression Eq. m4) can be represented in a form of a conventional Gaussian 



integral over the components of a complex 2A^— component vector S 

N 



Si 

S2 



Sp=l,2 — ( Si 



qiP) 



[dS] 



n 



dReSj dlmSj 



Det" 



Eb — flh 



Hi, 



[dSi][dS2] exp^ -S^ 



-i E, 



- Hin + iV 
kI 



J[dSi][dS2] exp |k(sIS2 - S|Si) + iEb{S\Si ~ S|S2) + z(Sl , S 

exp |iS7f,(SlSi + S|S2) - (Sl, S 




(66) 



At this point it is worth mentioning that all eigenvalues of Tie/ = -ffin — *r (scattering poles) must be situated in the 
lower half of the complex plane Imf < 0. Formally it is ensured by eigenvalues of the matrix F being real non-negative. 
We see that it is due to this fact that the Gaussian integral above is convergent | |103| (for O real other terms in the 
exponent are purely imaginary and do not spoil the convergency; at the end we may continue analytically). 



16 



The following comment is appropriate here. In principle, one can deal directly with the potential ^{E, Y, k), Eq.(|5£ 
and succeed in findi ng th e convergent Gaussian representation for the generating function everywhere in the complex 
plane E + iY{see [^,104|). However, the evaluation of the averaged generating function and subsequent restoration of 
the eigenvalue density turns out to be quite a dounting job. This is the reason why we decided to deal in our particular 
case with a less general, but more tractable representation Eq.(|66|) allowing to evaluate the generating function for 
two real parameters E, Q and then to continue analytically ^ —^ iY as explained above. 

Returning to our problem we represent the determinant in the numerator of the generating function Eq. (|6^) in the 



form of a Gaussian integral over a 27V— dimensional vector x — 



Xi 

X2 



whose elements x 



(p) 



j = l,2,...,7V;p=l,2 



are anticommuting (Grassmannian) variables (see reviews [100,101| for more details): 



Det" 



l)^^Det 



E - n - H„, - iT) { E + n - H,,n + ?T 

E + n- H„ 

kI 



kI 

-i(E-n- H,„ - if 



(67) 



(-1)^ J [dxi][dx2] exp {-«;(xlx2 + xlxi) + ix\ (e 



- H,n - iT] X2 + ixl{E + n- H,n + I 



.r)«} 



where [dx] — II^i ^xX^Xk- In contrast to the discussion above, the integration over Grassmann variables is always 
well defined and one does not encounter the convergency problem. 

Obviously, our generating function is the product of two Gaussian integrals defined in Eqs. (|66| , |67| ). It is convenient 
to introduce the notion of a supervector 



*2 



where 



1,2. 



Then one can write the generating function in the following " supersymmetric" form: 



where 



and 



Z{E,n-Eb,nb;n) = (-1)^^ / [d*] exp -(54*] + 5b/[*]) 



cS^*] = + i{n - r2b)*^/^fc* + i{E - Ei,)-^^Lkb^ + *M f ® Kl] * 



(68) 

(69) 
(70) 
(71) 



Before presenting the explicit expressions for the supermatrices L, A, Kb and E we would like to make a no- 
tational convention on arranging elements of supermatrices. All these (and subsequently appearing) supermatri- 
ces are assumed to act in the space of supervectors ^ whose element arrangement is defined in Eq.(|68|). Cor- 

Qii Qi2 



respondingly, we subdivide each 4x4 supermatrix Q into four blocks Q — 



Q21 Q2 



in such a way that 



^^Q^ = I]m,n=i *LOnm*n- cach of thcsc Q„„ blocks is in turn a 2 x 2 supermatrix Q,nn = 1 (mn) 



bb 



such that ^iQmn^n = St„Q{,r'S„ + SIqIJ''^^ + xLQ^r'Sn 



xlnQ^ff^^^ Xn- The indices 5,/ remind us of 



^m^bb 

"bosonic"/ "fermionic" nature of the commuting/ Grassmannian components of supervectors, respectively. 

It is necessary to note that in the present paper we use the same convention for Hermitian conjugation of 2 x 2 



supermatrices as in the paper [|101| : ( S.'^^ 



t 



This is different from the convention used in 



^bb ^fh 

-Qbf Q}f_ 

|100,105| and results in some subsequent differences in parametrizations. 
With these conventions the 4x4 supermatrices appearing in Eq.(^,^l]) are given by the following expressions: 



and S = 



I2 
k 



L = diag(l, 1, -1, 1); A ^ diag(l, 1, -1, -1) 
where I2 = diag(l, 1); k — diag(— 1, 1). 



i^h = diag(l, 0,1,0); 



(72) 
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B. Ensemble-averaged generating function. 



A quick inspection of Eqs.(^9|j7C|,|7^) makes it clear that the superintegral in hand is very similar to that emerging 
in the problem of ca lculation of the pair correlation function of eigenvalues of Hermitian random matrices at points 
E±il , see e.g. [101|. In fact, the two expressions coincide if one neglects the exponent Ss- The neglected exponent 5^ 
does not contain random variables and can not prevent us from using successfully the main steps of Efetov's standard 
procedure when evaluating the average value of the generating function. Below we give a short description of the 
main steps of the method; all further details can be found in the review [101| and references therein. 



• Ensemble averaging. One can easily perform the averaging over the Gaussian-distributed matrix elements of 
Hin by exploiting the identity 



(73) 



In order to write down the result of the ensemble averaging in a convenient form it is useful to introduce the 
supermatrix A with elements 



(nn) 
qq 



(74) 



where indices p and q are equal to b or / and we assumed the convention: = S> 



ensemble-averaged value of the corresponding exponent in eq. (p9[) can be written as 



exp 



Hi, 



L] * 



= cxp-^Stri2 



Now the 



(75) 



where the symbol Str stands for the graded trace StrQ — TrQbb ~ T^^^Qff- It is also useful to notice that 

an arbitary supermatrix Q. In particular, StrAA ~ ^^AL^. 

• Hubbard- Stratonovich transformation. 

As a result of ensemble averaging the superintegral representing the generating function ceased to be a Gaussian 
one. The further progress is based on the following identity: 



J [dR] exp I 



exp 



N 



^-Stri^ + mStrAA 

2N 



(76) 



[dR] exp { - — Stri?2 -|- iStiRA + 7VnStri?A| 



known as the Hubbard-Stratonovich (HS) identity. 

Now we can substitute this relation back into the averaged generating function, to change the order of integrations 
over the supervector ^ and the supermatrices R, and to calculate the corresponding (Gaussian) integral over 
exactly using the identity: 



j [d^] exp (-f^F*) = Sdef^i^, 



(77) 



where the notation Sdet stands for the graded determinant: SdetQ = exp Str In Q. It turns out, however, that 
in order to have both and J[dR\... convergent, one has to parametrize the set of supermatrices R in the 

following non-trivial fashion suggested in |100 9^ ( see detailed discussion in |101]): 



R = T-^PT- 



P = 



Pi - iSh 
P2+ iSi2 



Ptn Vrn 
^7n ^Qm 



(78) 



where Im, m = 1, 2 are real commuting variables, rj, rf are Grassmannians and the supermatrices T belong 
to the graded coset space U(l,l/2)/U(l/l)x U(l/1), and 6 is positive infinitesimal. 
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• Saddle-point calculation. 

Performing the program specified above one gets the following representation for the average generating function: 



/N " 
[dR] exp {-y Stri?2 + NnStrRA - Str In G} 

where 

G = Gi® In -iT® A; Gi = -IkT-l - {Eb - E)Kb + [n - 9.b)kbL -Eh-R 
and El — L^^^^t.L^^^^ . Now one can write: 



(79) 



(80) 



Str In G = TVStr In Gi + Str In 
The second term in this expression can be rewritten as: 



In - iT (g> (AG^^) 



M 



Strln 



/at - if ® (AG^') - ^ Str In {J4 - ^7a(AG^')} 



(81) 



which easily can be veryfied by expanding the logarithm into the series, exploiting the orthogonality condition, 
Eg. (p6[) , in each term of that expansion and resumming the whole series back |106 . Up to the present point we 
did not make use of any approximation and our calculation was essentially exact. However, we are particularly 
interested in the limiting case of many resonances iV ^ 1 coupled with few open channels M <C iV. In this 
limit we expect that the resonance widths are of the same order as the mean separation between adjacent 
resonances A cx 1/A'^. Therefore, we can restrict our attention to the case Y ~ 1/-/V, and, correspondingly, 
consider 1/N. The second fact that should be taken into account is that we are actually interested in the 
limit K, Eb — E , rib — ft when calculating the generating function. These facts taken together make it clear 
that it is sufficient to expand the logarithm Str In Gi in the exponent of Eq. (|7^) with respect to k, Eb — E,ilb — fl 
and retain (apart from the leading terms) only terms linear in these variables. At the same time we can just 
neglect all these variables in the term J^aLi Str In {/4 — i-faiAG^^} because of the condition M ^ N. As the 
result, we have: 



Z{E, n; Eb, nt; k)) = / [dR] exp -N£[R] + 5L 



where 



C[R] = -Stri?2 + str In {-EI - R) 
SC ^ NQStrRA + TVStr [^kSl + {Eb - E)Kb + {^b - ^)KbL] {-EI - R)-^ 



(82) 

(83) 
(84) 



The form of the integrand in Eq.(^) suggests that it can be effectively calculated by the saddle-point method 
exploiting the large parameter N 3> l.The saddle point equation is determined by stationarity of the "action" 
C[R\ and has the form R = {-EI - R)'^. At the same time the discussion above makes it clear that the terms 
entering SC are of the order of unity when N 00 and should be disregarded when seeking for the saddle-point 
solution. 

Actually, it turns out that there is a whole continuous manifold of the saddle point solutions Rs equally con- 
tributing to the integral Eq.( p2|) in the limit N oo:Rs = —E/2 + iiTVscT~^AT, where Vsc = i^sc{E) stands for 
the semicircular density, Eq. (||). Introducing a new set of matrices Q = —iT^^AT satisfying the conditions: 
= —I4; StrQ = we finally write down the averaged generating function as: 



{Z{E,n;Eb,nb;K)) ^ f [dQ] f[Sdet-' 

exp |-7V7ri^,cf^StrQ(A - KbL) - N-Kv.^VLb^trQKbL- (85) 
iirpNnStvQtL + N{Eb - E)[E + -KVscStvQKb]^ 
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To obtain the required function p{E^ - 



-if2)„ we should substitute this expression into the relation Eq.(f32l) 

d 



Upon 



doing this one immediately notices that the first term ^^'<^{EbMb)^{E,n) ^(■^) produces a contribution which is of 
the order of N, whereas the term ^ ^^^{Eb,nb)~^{E.n) M ^■^) produces a much larger contribution of the order of N^. 
Retaining only leading terms as long as ^ cxd we arrive at the following expression: 

p{E,-iu)e 1 



M 



i\laEk 



[dQ]Str(gA)Str Q(A - KbL) 



(86) 



ITTVsc^aQk 



exp 



-|strQA- |strgi]L 



where we introduced the "scaled" variables uj = 2'kVsc{E)NVL and e = 2tti'sc{E)Nk and the correspondingly rescaled 
p{E, -ifl)^ p{E, ~iuj),2T:Nvsc{E). 



C. Distribution of resonance widths: general expression. 

To perform the explicit evaluation of the superintegral on the right-hand side of Eq.(p6|) one has to employ the 
parametrization of the manifold of the Q— matrices suggested by Efetov [100|. To make the presentation self- 
contained we present such a parametrization in Appendix B. In the same Appendix we evalute also some supertraces, 
superdeterminants and combinations of the matrix elements entering Eq.(86) as well as other superintegrals we use 
later on. Upon substitution of these expressions into Eq. (^) one can perform the Grassmannian integration trivially 
and obtains: 

exp[iw(Ai - A2) + ie(| px | sin^i -Hi | ^2 | cos(/!i2)] TT ( , ,^ ) -F£,£;(Ai,2, 0i,2) 

where 

^'£,£(Ai,2,(/'i,2) = ie(Ai - A2)(| p\ I sin 01 -i\pi \ cos02) - ?;A2 x (88) 

{§(1 A^i I sin 01 -Hi I I cos 02) + [^(| P\ I sin0i ~i\p2 \ cos02)^ + 
HI Ml P - I Ai2 P -2i I A'iA^2 I sin (0i - 02))] } 

with g„ = 2/Ta - 1; A? = 1+ | p^ p, A^ = 1- | p^ p. 

Having in mind that actually we have to perform the analytic continuation — — > y in the right-hand side of the 
Eq.(^ we introduce two functions i^i,2(e, — «w) according to the following definitions: 

Fi{e,-ioj) 



— / dAi / d4>i exp [iujXi -I- ie | ^1 | sin0i] TT {ga + Xi)~^ = 
Ji Jo 

/ n ('^'^'^^^P-'^'^^'^) / dXiJo{€Jxl-l)exp[~Xi{~iuj + Y^Sa)] = 

•^0 a=l -^1 

poo 

/ '[[{dSaexp-Sa9a)<i>i{e,-iuj + ^Sa); $i(e,y) = 

•^0 a=l 



(89) 



and, similarly 



2 i-l I-2TT M 

F2{e, -iuj) = 1— / dX2 j d(f>2 exp [-iujX2 - e \ P2 \ cos 02] {ga + A2) 
^ J-i Jo 

M 1 

T.'^i^'^^a}^ / dA2A2^/o(e I P2 |)exp-zc.A2 



(90) 



fc=0 

M 



sinh y'e 
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where 



M 



M 



M 



a—1 a—1 h^a a,fc— 1 c^a.b 



(91) 



and Jo(z) and /o(^) = Jaiiz) stand for the Bessel functions. 

One can easily satisfy oneself that the right-hand side of Eq. ( ^ ) can be expressed in terms of the functions Fi (e, —iuj) 
and i^2(e, —iuj) and their derivatives. This gives us the possibiUty to perform the required analytic continuation easily 
and to restore the two-dimensional density p{E, y); y — 2ttVsc{E)NY , see Eq.(|6^) in the form: 



d_ f p{E,v) 
dy V i^sc{E) 



^ _ lini <! A ( 



[ dy 



d_F\ 
97 



d 



de 



d 



2de 



F, 



dFn 



de^ dy 



Fi 



de^ dy 



dFi d dF2 d^Fi dF2 d^ 9^2 
4:^ — ^ + -:r^-:r^ - Fi- 



de de dy dy"^ dy 



dy^ dy 



2Fi 



dy 
dF2 
dy 



(92) 



The limiting transition e ^ is performed with the help of the identities: 



d 



d^ 



- linie— $i(e,j/) = lini e-^ — $i(e, y) = 2S{y) 



(93) 



which, in turn, are consequences of the following representations for the S— function valid for an arbitrary integer k: 



S{y) = Ck lim 



(e2 +y2)(2fc+l)/2' 



Ci = 1/2; C2 = 4/3 etc. 



(94) 



It is useful also to note that if we substitute <i>2(e, y) for $i(e,y) in Eq.(p3|), this will produce zero instead of the 
(5-functions on the right-hand side. The same is true also for terms like limg^o e^^2'&i(e, y)- Using these observations 
we easily pick up all nonvanishing contributions to Eq.(p2|). Summing them up we get: 

dFi' 



dpiE^y)^d 



F^e- 



de 



dy Vsc{E) dy e-o 
that immediately allows us to restore the density p{E^ y) in the form: 

p{E,y) 



Vsc{E) 



POD ( M \ 

F2{e^Q,y) [] dSa ] S{y + Sa) exp - ^ g^S^ 

•^0 Va=l J a a 



(95) 



(96) 



It is clear, however, that for any positive y the 5- functional constraint in Eq.(g6|) is never satisfied and the right-hand 
side is identically zero. This result is of course just a consequence of the fact of absence of scattering poles in the 
upper half plane of complex energies. 

We therefore consider y < from now on, make the substitution y — s- —y to the previous equation and consider 
y > after that. The M-fold integral can be further simplified upon using the integral representation: d{u) = 
^ dk exp iku. Finally we arrive at the following expression: 



piE,y) 

Vsc{E) 



^i[{5},y]^2[{.9},y] 



where 



^ / n ^ ^ E n 

a—1 a—1 b^a 



9a - gb 



(97) 



(98) 



and 



M 



M 



-^2 [{9}, y\^^l ^ d\e-y' Ylig, + A) = J2{-l)'^vi"\g}^ 



sinh y 



k=0 



y 



(99) 
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and the functions '{.g} are defined in the Eq.(pl|). 

It is clear that the function pE{y) = ^^^(f:) the meaning of a distribution of (scaled) resonance widths for those 
resonances £j = Ej + lYj whose positions Ej fall into a narrow window 5E around the point E of the spectrum: 

^ Ne 

PE{y) = J^Y.^{y~ 2iTVsc{E)NY,). 
i=i 

Such a window should contain a lot of individual resonances: Ne ~ ^ ^ 1 in order to make the statistics represen- 
tative. On the other hand, it should be small in comparison with the total width of the spectrum (in our model given 
by the widths of the semicircle) to ensure that the local mean level spacing A(£^) is constant across this window. 

The expressions Eqs.( p7[]9^ ) provide us with the explicit formula for the local- in spectrum density of scattering 
poles for a generic open chaotic system with broken time-reversal invariance and constitutes the main result of the 
present section. 



D. Properties of the resonance width distribution. 

• Weak coupling versus strong coupling: from distribution to power-law behavior. 

When all scattering channels are equivalent, i.e. have equal transmission coefficients Ta = T (hence, equal 
9a — g) the distribution Eq.(p7|) can be represented in a quite simple and elegant form. Indeed, for this case we 
have: 

^ r°° 1 1 

where r(Af) = {M — 1)! stands for the Euler Gamma- function. We also can write the function !F2i9,y) in this 
case as 

^.(.. D-D- (1' (^) (101) 



k 

PE{y)can be written as 



where ( 1 stands for the binomial coefficient. One immediately sees that the expression for the density 



tl^yM-id^ ^„_,,sinhy 
r(M) ^ dy 



PE{y) - \^y''-'h^ [ e--= ) (102) 



Remembering, that the "weak coupling" limit corresponds to large values of the parameter ^ 1, one immedi- 
ately notices that the distribution of resonance widths is exponentially cut at ?/ ^ ymax = niaxajg"^}. Thus, in 
the weak coupling limit ymax ^ 1 and we can put ^ ~ 1 everywhere. This procedure immediately results in 
the well-known distribution for the case of equivalent channels, see Eq.( |102| ). The condition t/ <C 1 just means 
that the resonances are too narrow to overlap with each other. It is therefore natural that the distribution 
simply follows from a first order perturbation theory, see the discussion in the beginning of the present section. 

As long as the coupling to continuum becomes stronger, the parameters ga decrease towards unity. When one 
or more transmission coefficients Ta attain their maximal value Tq = 1 a drastic modification of the resonance 
width distribution occurs. Indeed, it is easy to see that when Ta (and hence the corresponding ga) tends to 
unity, the function [{g}, y] behaves proportionally to exp (— y) at large enough y. This decay exactly cancels 
the growing exponent exp (y) originating asymptotically from J-'2 [{g}, y]- As the result, the distribution function 
PE^y) must show a pure powerlaw decay in its tail, see Fig. 2. 

To determine this power explicitly one should make more accurate estimates of the asymptotic behavior of both 
J-i and J^2- Substituting sinhy w ^e^ in the definition of J^2 one can write: 

:F2[{9},y^^] = leyJ2t^U„ U, = J2{-lf-^V^^'^ (103) 

p=0 ^ fe=0 ^ 
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Using the definitions of the coefficients 2?^*^'' , see Eq. (^) one finds that 

M M M 

f/o=n(5a-i); c/i=-En(5fc-i); - u,^{~irp\ n (^c-i) 

a— 1 a—lb^a ai ,. . . .a/ — 1 c^^ai . . .a; 

We see that the leading power of y in the asymptotic behavior of the function J-2 [{g}-, y] is essentiaUy determined 
by the niimber of parameters ga which are simultaneously equal to unity. Let us suppose, for definiteness, that 
exactly / quantities 51,(72, ■■■,gi are equal to unity, whereas all other M — I parameters gi<a<M are larger (and, 
for simplicity, are all different). Then Uq = Ui — ... — Ui-i = and the leading behavior of the function 
^2 [{g},y] is given by: 

y 

J-2[{5},2/-oo] = (-l)'n2^ \{{ga-l) (104) 

^ a=l + l 

Under the same conditions one can determine the leading asymptotic behavior of the function J-i [{g},y\ by 
calculating the integral in eq.(p8|). The contribution of the fold pole aX k = —i gives: 

[{.,}, . ^ 00] . ^y^-^e-y ^ ^ (105) 



This finally results in the desirable asymptotic decay law for the distribution of resonance widths: 

I 

2y 



PEiy^(x,) = — (106) 



The physical origin of such a tail in the width distribution is discussed below. It is interesting to note that such 
a behavior means that the positive moments of the width distribution dyy^ pE{y) are apparently divergent 
as long as A; > 1. 

It is also instructive to consider briefiy the particular case of very many equivale nt ch annels: M ^ 1; ga ^ g 
for any channel. We find that it is most convenient to rewrite the expression Eq.(102) in an equivalent form: 



1 



2r(A/)2/2 



PEiy) ^ oT^^.n 2 / dtexp-(t-Aflnt) (107) 



and to evaluate the integral by the saddle-point method. The exponent is maximal in the vicinity of ts ~ M . 
When this point is inside the integration region, i.e. y{g~ 1) < ts < y{g + V) we have a nonvanishing contribution 
to the integral. In the opposite situation the density of resonances vanishes exponentially when M ^ 1. Picking 
up the nonvanishing contribution we obtain: 



M _M_ <- „ JLL 

J SP'' g+l ^ y ^ g-1 



Two conclusions can be drawn from this expression: 

1. In the limit of large number of channels M ^ 1 the distribution of resonance widths shows a gap: there 
are no resonances with widths smaller than y„i = 

2. A region of power-law behavior psiy) M/y"^ exists not only for the critical coupling g = \, but also in 
the vicinity of the critical point: g — 1 <C g. However, only for (7=1 the power-law domain extends to 
infinity. 

The formation of a gap (a strip in the complex energy plane free of resonances) was first noticed in the numerical 
experiments by Moldauer ]107| ] long ago and discussed in much detail by Sokolov and Zelevinsky Q later on. 
Gaspard and collaborators (see references in [|j) observed such a gap in their studies of chaotic scattering in 
the so-called three-disk systems. Semiclassically, the number of open channels was very large and comparable 
with the number of resonances. In the limit M, N 00 but m — M/N finite the expression for the resonance 



width distribution was obtained and analyzed in the papers l5qJ5q|. Our expression Eq.(lOS) obtained under 



the conditions 1 << M << iV perfectly matches the m << 1 limiting case of their expression. 
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• Mean resonance width: the Moldauer-Simonius relation. 

Having at our disposal the explicit formula Eg. (p7|) we can, in particular, easily calculate the mean resonance 
width: 



(y) 



(109) 



pOQ -X poo -I poo 

/ dyyTiT2 = -— dXT[{ga + X) dkY[ -/ dyye-y 

Remembering the relation between ga and the transmission coefficients Tq, see Eg. (|50|) , and using the definition 
of the scaled level width y = ttF/A we can represent the last result in the form of a relation between the mean 
resonance widths (F) and the transmission coefficients Ta'. 

A 

{r)=--J2\n{l-T,) (110) 



2tt 

a=l 



which can be also rewritten as: 



M 



I lliSa) hcxp-^ (111) 

a=l 

in view of the definition Ta = 1— | (Sa) 

The latter formula is well-known for a long time in nuclear physics as the Moldauer-Simonius relation [108|. 
It was derived for systems with unbroken time-reversal symmetry by averaging the 5— matrix over the energy 
spectrum and using the unitarity condition. The fact that we recovered this relat ion by ensemble averaging is 



in good agreement with the well-known ergodicity of the Gaussian ensembles [109|. The logarithmic divergency 
of (r) a t the critical coupling Ta = 1 is a direct conseguence of the 1 /y^ decrease of the probability distribution, 



see Eg. (11061). 



• Strong chaos as the origin of the power-law tail of the resonance width distribution. The results 
presented above suggest that the powerlaw decrease 1/y^ should be typical for chaotic systems strongly coupled 
to continua and is one of the clear manifestations of the strong overlap between individual resonances. It is 
therefore natural to try to understand the origin of such a tail gualitatively in terms of the underlying chaotic 
dynamics. 

To this end it is interesting to mention that a little different powerlaw distribution of resonance widths, that 
of the form p{y) oc y^^^^, was observed in numerical studies of a guantum chaotic system with (guasi) one- 



dimensional "diffusive" dynamics in the case of strong coupling to continua [110|. The authors suggested a 
transparent gualitative explanation of this effect based on the fact that the resonance width is proportional to 
the inverse life-time for a wavepacket injected into the system. The latter is determined by the time of the 
classical diffusion: tdif = D/L^, where L is the distance from a given point to the closest (strongly absorbing) 
boundary and Z? is a classical diffusion constant. For a semi- infinite sample this reasoning immediately gives a 
powerlaw width distribution p(F) cx ^ c>c F^"^/^. 

Let us remind that the present model is based on the use of the Gaussian random matrices. Physically, it 
corresponds to the case of strongly chaotic classical dynamics for the closed system |2^,Q . For such systems 
there is a typical time scale 6t determined by an inverse Lyapunov exponent after which the system effectively 
loses a memory about its initial conditions and can be found in any part of the available phase space on the 
energy shell with egual probability. Such systems are known as the ergodic ones. Let us show that it is just 
that type of classical dynamics which is responsible for the powerlaw tail l/y^ of the width distribution. 

To understand this fact let us consider as a particular, but generic example: a particle moving with a velocity 
V inside an irregular-shaped two-dimensional cavity of area A and circumference C oc A^^^. The chaoticity is 
considered to be so strong that complete " loss of memory" occurs after few collisions with walls so that St can 
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be estimated as St cx C/v oc A^l"^ jv. Let us make a small opening of the width d « C in the walls so that the 
particle can escape from the cavity whenever it hits the opening. Subdividing the observation time in intervals of 
the order of 5t we conclude that the probability to escape during one interval is just po = d/C cx c?A^^/^ << 1 in 
view of the ergodicity and escape events during the subsequent intervals can be treated as independent (memory 
loss). Then the probability to stay inside the cavity for a large time and then to leave within the interval 
[ie, + dt] can be estimated as P{te)dt = exp —pot^/St. 

Considering our system semiclassically, we associate the wavelength \d = with our particle and can estimate 
the number of the quantum mechanical states available inside the closed cavity as iV cx Since the energy 

E = the corresponding mean level spacing being of the order of A = E/N is proportional to mv'^X^/A. 

Interpreting the inverse escape time h,t~^ as the resonance width F and measuring it in units of A: y = F/A, 
one can find the distribution of y to be given by: P{y) = ^^exp—^^^ w for y ^ Msc, where Msc = 

^ mv^X^ A^i^ — d/\d: which coincides with the quasiclassical estimate for the number of open 
channels for the present problem. 

We conclude that the semiclassical arguments faithfully reproduce the same powerlaw tail of the resonance width 
distribution as that obtained from our random matrix model. Therefore, we expect such a tail to be a universal 
characteristic of the chaotic quantum scattering problem independent on the specific details of the underlying 
classical dynamics being sufficiently chaotic to ensure an exponential escape from the compact scattering region. 
The following comment is appropriate here. Dealing with realistic models of open chaotic systems containing 
no random parameters one always performs statistics over an interval of energies 5E on a real axis containing 
many resonances: 5E ^ A, but being small enough for a systematic variation of the smoothed level density 
y{E) to be neglected: 5E ^vj^ (c.f definition of the quantity pE{y))- One may expect that universal 
features of such statistics are adequately reproduced within the framework of the stochastic approach, but only 
on the level of " local-in-spectrum" characteristics calculated at fixed value of E. Indeed, any spectral averaging 
in the stochastic model performed on a scale comparable with the radius of the semicircle unavoidingly mixes 
up data corresponding to very different values of the transmission coefficients, the procedure washing out any 
relevant physical information. In particular, it seems quite meaningless to consider quantities like the "globally" 
averaged resonance width Tgi — jj X^feLi Tfc, where the summation goes over all N resonances. In our model 
this quantity can be trivially found from the sum rule: NTgi — — 2Tr luiH^f — 2 7^. and can not be related 
to any particular transmission coefficient. This fact, however, should not be misinterpreted as impossibility to 
have universal statistics of S*— matrix poles within the stochastic approach as discussed in pa]. Rather, the 



quantity F^; can be found via the direct integration of the universal local expression (F) , Eq.(llC) over the 
energy E, upon substituting ther e th e energydependent values Ta{E), A{E) from Eq. (^9|) . Indeed, the following 
integral can be easily evaluated [ 1 1 1 1 : 



/•2 

lim y^Tk^N / {T[E))u[E)dE = (112) 



N 

k 

M 



resulting in the expected expression 2 7a as long as all 7a < 1. 

We see that the result of the integration is always finite for any 7^, thus concealing a specific role of the critical 
coupling 7a = 1 when resonances with divergent local mean width occur sufficiently close to the center of the 
spectrum. 

Of course, taken literally this divergency has sense only in the limit of infinite number of resonances iV — > cx3. 
For any finite N all resonance widths are finite and in any ca se c an not exceed the upper bound 2^^ 7a. 
Basically, it is related to the fact that the distributions Eqs.( |97 , [l02 ) cease to be valid for the domain of very 



broad resonances having widths F ~ 1 (correspondingly, y ~ N). Alternatively, one may say that for large, but 
finite N the Moldauer-Simonius relation is to be modified in a narrow domain 5E cx 1/A'^ in the vicinity of the 
energy E — Q, see M. 



Expression Eq, 112 can be also used to describe an interesting phenomenon happening when some coupling 



constants 7a ( e.g. for the channels a = 1,2,..., Mi < M) exceed the critical value 7 = 1. Under this 
condition the result of integration in Eq. (|ll2| ) is less than the exact sum rule value 2 X^a^ la by the quan- 
tity 6T = 2^^^^ (7a — 7^^) > 0. This deficit reflects the existence of Mi "broad resonances" of the widths 
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Ta = 2 (ja — ^~^) » A cx 1/N; a = l,...,Mi. Such resonances can not be described by the distributions 
Eqs.(97 102 ) which cover only the resonances satisfying y = ttF/A < oo in the limit N oo. When coupling 
constants ja increase further the group of broad resonances moves away from the real axis, accumulating in the 
limit 7a ^ 1 the lion's share of the total widths Tgi. The remaining N — Mi resonances become progressively 
more and more narrow and their widths arc well described by the distribution Eq.(^). 

The effect of reorganization of the resonances into two essentially different groups due to strong coupling to 
continua ( the "trapping phenomenon" ) was studied in some details in the papers ]53|,|55|,p9|,p0[ . 



IV. STATISTICS OF SCATTERING PHASE SHIFTS AND TIME DELAYS 



In the present section we are going to study in much detail the statistics of individual scattering phase shifts 9a 
and their derivatives, both over the energy E and over an arbitrary external parameter X. As we already mentioned 
in Sec. 11 (see Eg.^), it is convenient to characterize the phase shift properties via the spectral density: pe.x{z) = 

17 X^ali '^(■^ ~ Za{E,X)) of eigenvalues of the matrix Kx{E) = ttW^ (^E — Hin{X)j W. Our object of primary 
interest is the correlation function 

K.EM,x{zi, 22) = {PE-n/2,~X/2{zi)pE+n/2,X/2{z2)) - {PE-n/2-X/2{zi)){pE+n/2,X/2iz2)) (113) 

knowledge of which, in particular, allows one to study statistics of "partial delay times" Ta = and also the 
corresponding parametric derivatives. 

Before addressing the issue of the spectral density correlations, it is instructive to consider in some detail the 
calculation of the average sp ectr al density {pe.q{z)) for the few channel case. This quantity is less informative than 
the correlation function Eq.( |ll3|) , but that simple calculation serves as a reference point for more interesting cases. 
Let us mention, that in the limit M oc 3> 1 the phase shift density was found earlier by Lchmann and Sommers 



A. Averaged spectral density of K— matrix 

The averaged density can be easily found provided the following functions are known 

ftxi^) = U—-^^-^) (114) 
\ z±ie- Kx[E) / 

in view of the obvious relation: pe.x{z) = lim^^o Ini/^ restrict our attention in the present context by 

Ie x=oi-^)' omitting all the indicesi, i?, X for the sake of brevity. 
The function f(z) can be formally written as: 

n Det ({z + J)iM - k) 
/(z)= lim— (InZ(J)); Z{J) = ^^-^ ^ (115) 



Here and below we imply z = z — ie for the sake of brevity, implying f{z) to be analytic in the lower z half-plane. 

Due to the normalisation condition Z(J = 0) = 1 one can write f{z) = limj^o '§j{'^{J))- order to perform the 
ensemble average of the generating function Z{J) in a standard way one should first get rid of the following unpleasant 
feature: the random matrix Hin enters the expression for the generating function only via the matrix K. To this end, 
we can use the identity Eq.(^) and write down the determinant in the denominator of the generating function as: 



z 



*^Det [Im - z-^t:W\E - Hiny^W^ = z*^Dct [In - z-^7r{E - H„,y^WW^j = (116) 

z^'Dct-^ (e- Hin) Det - i?,„ - ttz'^WwA 



After performing a similar manipulation with the numerator of the generating function we can write: 
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z + J 



M 



T)et[E - Heff{z + J)] 



(117) 



where we introduced the notation Hf,ff{z) = Hin + jWW^ . 

Now we can use a standard procedure and represent the determinants in the denominator/numerator of the preceed- 
ing equation by Gaussian integrals over N commuting/anticommuting variables. After introducing 2iV— component 
S, 

X 



supervector 'S' = 



I , we have: 

Dct[E~Heff{z + J)] 
Dct[E~H,ff(z)] 



[d^] exp |-i\I>^ 



{E — Hi, 



T(E)U 



(118) 



where the supermatrix U = diag(z~^, (z + J)^^) and f = nWW^ , as before. 

Now one trivially performs the averaging over the ensemble, see eq.(73) and decouples the emerging "quartic term" 
in the exponent with help of the Hubbard-Stratonovich transformation. In the present simple case such a decoupling 
is possible when one uses the set of 2 x 2 matrices P defined like in Eq.(^^ as the integration manifold. Changing 
the order of integrations, performing the integration over explicitly, and copying steps used to derive Eq.(|8l]), one 
obtains: 



/ Det[E-H,ffiz + J)] 
\ Bet[E-H,ff{z)] 



M 



dPe-^^(P^ J] Sdet-i h - jaU{E - P)-' 



(119) 



where 



C{P) = 2 ^^''^^ + ^^"^ " ^) 



Now it is convenient to use that: 



Sdet^^t/. When combined with the preceeding expression it gives: 

M 



J{z) = lim ^ / dPe-^^(^) n Sdct-^ [iJ-' la{E Py' 



(120) 



The integral over P is calculated in the limit ^ M by the saddle-point method, with unique diagonal saddle- 
pointPs = {E/2 + i'KVsc)l2 accessible by allowed contour deformation for | |< 2. This immediately yields: 



M 



/(^) = lini A J]^ + - la{E/2 + invsc) 



J^o dJ 



a=l 



z - 7a(i?/2 + iirvsc) 



M 

V t 

^ z-7a(£;/2 + i7ri/sc) 



(121) 



analytic in the lower half-plane, from where we find that the mean density of eigenvalues for the matrix K is given 
by a sum of Lorentzians (with z real): 



1 ^ 



(122) 



For the particular case of one open channel M = 1 the Lorentzian form of the average spectral densiy was first found 
by Mello . Actually, in that particular case one can check the expression Eq. (^8|) for the averaged 5- matrix using 
Eq.(122). Indeed, for M = 1 the 5*— matrix is reduced to the only number S — exp (— 2iarctanz) = 1 — 2iz/(l + iz). 
Therefore 



(5) 



l-2i 



dzpE{z) 



1 + iz 



(123) 



The integrand has the only pole z_ = E^/2 — invscj in the lower half plane Imz < 0, and the corresponding residue 
immediately gives: 



(S") = 1 - 2i{-2TTi) 



Z_ - E-f/2 - ITlVscl 1 



1- i(^£; + V4-S2) 
1-f i(zP + V4-S2) 



(124) 



in complete agreement with Eq.(Es). 
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V. PAIR CORRELATION FUNCTION OF A'-MATRIX SPECTRAL DENSITIES 



Let us now turn our attention to the calculation of the pair correlation function K,E,n,x{zi, Z2), Eq.(113). To this 
end let us introduce the function 



f{zi,Z2) 



Tr- 



1 



-Tr- 



1 



zi-ie- K_x/2{E ~ n/2) Z2 + ie- Kx,2{E + 



(125) 



related to the correlation function in Eq.(113) as: 

1 



E.n.x 



{ZI,Z2) 



27r2M2 



Re/c(2:i, Z2); fc{zi,Z2) = f{zi,Z2) - f (Zi)/+(Z2) 



(126) 



Pe rforming with each of the two traces of resolvents in the Eq.(125) the same manipulations as presented in 
Eqs. (115-117) one obviously obtains the following representation: 



f{zi,Z2) = 



32 



dJidJ2 



7(1) 7(2) 



^ ^ _ / Det[i; - n/2 - ge//(-X/2; zf)]Dei[E + r>/2 -~ H,ff{X/2; Z)'')] 
\Det[i? - n/2 - H^ff i-X/2; Zy^,)]Det[E + n/2 - H,ff{X/2; zf^,)] , 



M 



^(Jl, J2) 



|Jl = J2=0 



where we introduced the notations: Z^^ ^ Zp + i(— 1)^6 + Jp; p = 1,2 and H(.f f{X; Z^f) — Hm r -^Tpi-^-^in 
-^WW^. 

This expression is quite close in its form to the generating function Eq(|^) appearing in the calculation of resonance 
widths distributions and we can use a similar representation for it in terms of the Gaussian(super) integrals (cf.Eq.dss" 



(127) 



:F{Ji,J2) = i-'^ f J M*] exp i^-iE-if^ ~ LA'i/ + i-i'^t (g) (lC/) 



(128) 



exp i i*^ ( Hin L ) - i 



X 



'2A^i/2 



L] * 



where IJ-^ = diag(Z^^2o, zf\zfl^^, zf) = ^^(/^ + A) + f (/4 - A) + diag(0, Ji, 0, J2) and notations for the super- 
matrices L, A and the supervector ^ are the same as in Eqs.( |68| , [72|) 



The subsequent procedure of dealing with the ensemble average in Eq.(12S ) is exactly the same as that presented 
in details in Sec. III. The only difference is that the average is performed only over the GUE matrix whereas the 
matrix is considered to be arbitrary, but fixed from the same ensemble, see [p8|. As a result, one has (cf.Eq.^9|): 



T{,h,J2)= / [di?]exp<j-^ 



Stri?^ - f7Stri?A 



- StrlnG-^ 



(129) 



where 



so that 



G-^ = Gf - r ® [/; Gf = 



{Eh -R)(E)l 



27V1/2 V " 



N 



.A 



(130) 



Str In = Str In Gf + Str In j/jv - (e «) C/) 
In turn, one can expand Gf^ in a series with respect to X : 



Str In Gf = A^Str In [Ei^ - -Rj " XI ^ ~ 

1=1 



I 



H, 



(1) 



Gi 



Str 



AiEIi - i?)~i 



(131) 



(132) 
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Now we use the fact that for any typical GUE matrix holds: Tr 



H 



(1) 



2p 



= 0(iV); Tr 



H, 



(1) 



2p+l 



= 0(1); where p > 



is an integer. It is therefore evident, that in the limit N oo the only nonvanishing term in the expansion above is 



that with 1 — 2. We also can put effectively Gf ~ yEI^ -~ Rj 'S) In in the second term in Eg. ( 131 ) and represent it 

in a form of a sum over channels, see Eq.(|8l|). 
Collecting all the relevant terms, we obtain: 



exp 



^(Ji, J2) 



[dR] exp 



-Str 



N 

M 



-Str7?2 + Strln {Eh - R) 



(133) 



A(£;/4 - Ry^ [ Yl Sdet^^ h - laU{Eii - Ry^ 



This integral can be evaluated in the usual manner by saddle point method in the limit N ^ 1. The saddle 

-1 



point manifold is parametrized again by the matrices R = — ttv^cQ = 

z(i) 7(2) 

Sdet?7 



EL-R 



Remembering also that 



7(1) 7(2) 

^ J=0^ J=0 

integral over the graded coset space 



, we obtain the following representation for the correlation function /(zi,Z2) in terms of the 



lim 



92 



Ji^0,,72^0 dJidJ2 



M 



[dQ] W Sdet"^ 



a=l 



U^^ - 7a ( —ii - TTVscQ 



(134) 



exp <j -^StrQA + ^StrQAgAl 
2 8 I 



where we introduced scaled variables: uj = TrVscNfl; x — irVscX. 

Remembering the definition of the supermatrix U^^ and performing the expansion of the superdeterminants up to 
the second order with respect to Ji, J2 one finds: 



lim 



92 



M 



dJidJ2 



7a 



E 



-I4 - TTVs 



(135) 



where 



Ba = ^(/4 + A) + ^{h - A) - 7a ( - TTi^scQ ) ; Ci= diag(0, 1, 0, 0); C2 = diag(0, 0, 0, 1 



^2 , 



E 



(136) 



In order to evalu ate the in tegral over the coset space explicitly we substitute the corresponding expressions (see 
Appendix B, Eqs.(B12-B15)) in Eq.(135) and perform the Grassmannian integrations remembering that in the chosen 
parametrization , see Appendix B, a nonvanishing contribution comes (apart from the terms proportional to the 
combination a*aP*P) also from terms in th e integr and containing no Grasmannians at all (the so-called Parisi- 
Sourlas-Efetov-Wegner (PSEW)theorem, see |105,101|). We therefore find: 



M 



M 

E 



1 



Z2(a) + i% ^ zi{b) - i% 



d\id\2 
-I (Ai-Aa)- 



■ exp 



-^^[X,-\2)--{\\-\l) 



M 

n 



-DM 



a—1 ^ ' 



(137) 



■ M 

E 



22(a) + i7aAi 2:2(0) + i7aA2 



2?/ (a) 



M 

E 

b=l 



21 (6) - «7bAi 21 (6) - i7bA2 



2?6(&) 



M 

E 



Vlia) ^ V}{a) 



where the notations I?/^b(a), 7q, 21^2(0) are explained in Appendix B. 
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It is easy to see that the expression in the first hne of Eq.(137) is just the so-cahed "disconnected" part / (zi)/+(z2) 
of the corresponding correlation function, which is given by: 



y i y I 

^ 2:2(0) + 1% f-[ zi{b) ~ ijb 



= ( Tr- 



1 



Tr- 



1 



Z2+ie-Kx/2{E + n/2) 



The "connected" part of the correlation function fc{zi, Z2) = f{zi, Z2) ~ f {zi)f^{z2) can be written in the most 
elegant form by noticing that: 



Z2{a) + i%Xi Z2{a)+i%X2 



and 



I?b(a) 

zi{a) - 

ll I A^i ? 



Zi{a) - i^a\2 
ll I P _ 



dzi Vbia) 
dz2 Vbia) 



(138) 



Vl(a 



In- 



(a) 'D^ia) dzidz2 Vb{a) 



Taking these relations into account one finally obtains the following compact expression: 



fcizi,Z2) = 



^ dXidX2 



-^LuiX,-X2)--i\l-Xl) 



(139) 



(140) 



dzidz2 



M r 

n 



zi{a)z2{a) + i^aX2izi{a) - 02(a)) + 7a 



zi{a)z2ia) + i%Xi{zi{a) - 22(0)) + 7^ 



with 7rz^sc7a and Zp{a) — Zp — jaE/2. This expression constitutes one of the central results of the present paper. 
In the rest of the present section we are going to use this relation intensively for extracting statistical properties of 
scattering phase shifts and their derivatives. 



A. Correlations of phase shift densities at fixed values of energy E and external parameter X 



The general expression eq.(140) can be further simplified in the particular case lu ^ x = 0. Physically this means 
that we are interested in studying correlation of phase shift densities at two points 61 and 62, but at fixed values of the 
energy E and the external parameter X . Let us further assume that all channels are equivalent 7a = 7; a = 1,2, M 
for the sake of simplicity. Introducing notations: A = z\Z2 + 7^; -B = 7a (z\ — 22) where zi.2 = 21.2 — -E'7/2 we can 
write: 



/; 



U—X — ^ 



{Z\,Z2) 



dXi- 



92 



1 



dX2 



dzidz2 {A + iBXiY^ J -I (Ai - A2)2 
The integration over A2 can be easily performed yielding: 



{A + iBX2) 



M 



/, 



LO — X — Q 



{ZI,Z2) 
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M 

■E 



M \ i-iB) 



dz,dz2^^ \ I J {I -I) J, [A + iBXiY 



dXi 
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dzidz2 



M 

E 



M\ (-1)' //-1\ (-1)« 



1) ^ 

g=0 



{(Ai + l)'-^-(Ai-l)'-i} 
A-iB'^'' 



(141) 



(142) 



A + iB 



- 1 



where in the second line we expanded the brackets (Ai ± 1)' and performed the remaining integration over Ai explicitly 
in each term. 

After differentiation over zi , Z2 with help of the relations 



^ dB ^ dA 
A- B- 



dzi,2 



dzi 



7(7^ + z2); A' + B'^iY + ~zf)iY + Si) 
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the sums over /, q can be performed explicitly as well. As the result, we obtain: 

/r^=°(zi,z2) 





'A-iB' 




( 








A + iB 



(143) 



In this point it is convenient to pass from the variables zi , Z2 to new " angular" variables 9i , 02 defined as: 

2i = 7 tan 6*1, Z2 = 7 tan 6*2 (144) 

We obviously have: 

,, sin(g'i-g2) A-tB ~ ~ 

B^T ~ —'^ A , = exp~2z(6>i ~ 6*2) 



cos Ui cos t/2 



Remembering the relation Eq.(126) between the spectral correlation function ICE,n,xizi, Z2) and /(zi,Z2)we notice 
that the pair spectral correlation function of the densities of the "angles" 6 defined as: 



ic{0i,e2) = {pe{0i)pe{02)) - {PE{ei)){PE{02)) 



(145) 



1 ( 

(^) — ^ (5 [ fii - arctan 



PeK 

can be written in a very simple form: 



E/2 



sinM(0"i - §2) 
TrMsin {§1 - 02) 



(146) 



One immediately recognizes in Eq.(146) the pair correlation function of the Dyson Circular Unitary Ensemble ( see 
the book where this object is called "two- level cluster function"). It corresponds to the following joint probability 

density of M variables 0a = arctan - E/2\ a = 1, 2, M: 

Pm(6'i, ...,0m) ^ const x J| 



^2i8a _ ^2i6^ |2^ 



e " — e ° I cx 

a<5 a<b 



(147) 



with — 7r/2 < 0a < 7r/2. Assuming this probability density being proven, the joint probability density of phase shifts 
0a related to "angles" 0a as tan6'a/2 — —7 ^tan^a + , where E = E/{'kVsc) is given by: 



M 

'Pm{0i, 6**/) = 'Pm{0i, ^m) JJ^ 

a=l 



d0a 



d0n 



d0a 



1 COS^ 0a 

2^cos2 6*0/2 



(148) 



On the other hand one can write: 



sin^ {0a - 0b) ^ cos^ 0a cos^ ^{,(tan - tan 0i,) 



cos 



cos^ f^b 



7C0S2 6'a/2 7COs2 6*6/2 



sin^ (0J2 - 0t/2) 



M 



SO that using the identity: ria<fc ^aUb = (Yia '^a)™ ^ one obtains: 

Pm(^i, ^m) « i n (^'^/2 - ^^/2) fn ] (149) 

Using the relation between 0c and 0c and definitions of the quantities E, 7 one finds after a simple algebra that: 
cos2 0r 2 



C0S2 6lc/2 (l ^ il ^ J:^ 



2i7 ' 2 



4 7^ 



+ e- 



2i7 ' 2 



4 7^ 
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(7^(4 - E'')/ [1 + 7' + 7V4 - E-^\ ) 
(1- (5)*e»e) (1- (5)e-»e) 



where we made use of Eq.(124). Thus, we arrive finally at the following expression: 



M 



-2M 



(150) 



a<b 



c=l 



which is nothing other but the Poisson's kernel distribution, Eq.(|j). Here the phase shifts 9a m ay b e restricted to an 
interval < 6a < 2tt. Inve rting the argumentation, we prove, that our correlation function, Eq.(146) follows from the 
Poisson's kernel, Eq.( 150); this we have shown in the case of equivalent channels. 



B. Distribution of partial delay times and parametric derivatives of phase shifts 



The knowledge of the spectral correlation fimction Eq.(113) allows one to determine the distribution P^(T)of partial 
delay times 



dOaiE) 



dza{E) 



a = 1, ...,M 



dE l + Za{Ef dE 

The distribution Vt{t) can be easily found if one knows the joint probability density Ve{z, v) defined as: 

VE{z,v)^^(f:S{z~Za)s(v^^'^^''^^ ^ 



because of the relation: 



00 /"OO 



dE 



dzdyVsiz, v)d r + 



2v 



1 + 



(151) 



(152) 



(153) 



where angular brackets stand for the ensemble average as before [113|. 

To determine the joint p roba bility density Ve{z,v) we use its relation to the spectral correlation function 
ICE.n,x{zi, Z2) defined in Eq. (|ll3| ): 

'Pe{z,v) = M Jim.^niCE.n^x=o{zi = z- vn/2,Z2 = z + vn/2) (154) 



relations of this kind were first used in [114| and later on in 1 115 -117]. To understand its origin we write for small 
positive Q: 



E.n.x. 



=o{zi — z — v^/2, Z2 = z + vCl/2) 



M2 

M 



(155) 



M2 



,6=1 



Expanding Za{E - n/2) - Zb{E + n/2) at small n as Za{E) - Zb{E) - f + + ■•■ we immediately see 

that in the limit — > +0 a nonvanishing contribution toEq.( 155) comes from the terms with equal channel indices 
a — b and the resulting expression is equivalent to Eq.(|153|). 



To perform the limit +0 in the most economic manner we pass from the variables zi,Z2 to z = {zi + 

Z2)/2; Vs = {z2 — zi)/2uj, with Vg^w being the scaled variables Vg — wA/27r; uj = TrQ,/A. Correspondingly, the 
correlation function fc{zi,Z2) acquires the form: 
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fiz,Vs 



/OO pi 
dXi / dX2 Gxp —iuj(Xi 



A2) 



1 

(Al-A2)2 



(156) 



M r 



n 

a=l L 



- ^aEz + 7a - 2i'KVsc^aijj\lVs - v'ju}'^ 



According to the general discussion presented above we should take the real part of this expression and look for 
the term proportional to l/cj at w — > +0. It is easy to understand that such a singularity comes from that part of 
the integration region over the "non-compact" variable Ai, where Ai oc ui~^. After a natural rescaling Ai — t/uj one 
can easily extract the corresponding singular term, which turns out to come from the term with the second derivative 
in the expression above. Performing the calculation explicitly, we find the expression for the joint probability 
density of variables z and Vs 



Ve{z.,Vs) = - 



(157) 



where Va = [z^ — jaEz + 7^] /2'Ki'sc7a and we made use of the identity: Re dtf{it) = i dtf{it). 
Performing the integration in the expression above (cf. Eq.(B8[)) one obtains: 



VEiz^Vs 



(-1) 



M-l 



M 

n 



laEz 



-Vs) X 



(158) 



M 

^exp- 
6=1 



z^ - -ibEz + 7^^ 



21 



2TTI^sc7bVs 



n 



z^ - jbEz + 



icEz + rc 



-I -1 



2'KVsclb 



where Q{x) = 1 if x > and Q{x) = otherwise. We took into account that z^ — "faEz + 7^ > as long as | |< 2, 
which is just the case we are interested in ( we remind that the semicircle density VsciE) is non- vanishing for | E |< 2). 
Substituting the expression Eq.(158) into Eq.(|153|) one trivially performs the integration over v because of the S— 



function and obtains for the distribution function of scaled partial delay times Ts = ^ the following expression: 



-PAts) 



(-1) 



M-l 



27rMT? 



2dz 



M r 



(-1) 



il + z' 



M-l 



n 



laEz + 72 



M 



^exp- 



b=l 



z^ ~ -^bEz + -fl 



7ri^sc76Ts(l + z"^) 



n 

c^b 



jbEz + 7^ z^ ~ -fcEz + 7^ 



2ttMt3 



d<P 



M 



27ri/sc7b(l + 27ri'sc7c(l + 2^2) 

RbW 



(159) 



a=l 



M 

^exp- 

b=l 



l[[Rbi<j>) - Rem- 

c^b 



Here 



Ra{(f') 



(1 +7a) - laE silKf) + {jI - l)cOS0 



2TTI^sc7a 

and we changed the i ntegr ation variable: z — tan {(j)/2). 

The expression Eq. ( [159| ) provides the distribution of partial delay times for a general case of non-equivalent channels 
and thus constitutes one of the most important results of the present subsection. 

Further simplifications are possible if we restrict our attention to the particular case of equivalent channels 7^ = 
7; a= l,2,...,Af: 

M 



1 



2ttM\t"-^ 



(1+7^) - ^Esm(j)+ (7^ - l)cos( 



1 



exp ' 



2'KVscl 

[(1 -I- 7^) - 7£;sin0 + (7^ - l)cos( 



(160) 



2'n:VsclTi 

The last expression can be put in a more elegant form upon using the identity 



f{pcos(j) + qsin 



= 2 



cos 
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and introducing the quantity g = (7+7 
so that ^^/W^WThW = ^fg' 



'^^j^-KVsc (related to the transmission coefficient as g 
— 1. This finally gives: 



: 2T-1-1, see Eq.(|o|) 



^r(T.) 



M 



1 cos 

(_1)M 



exp 



1 



A/ 



g + \J g"^ — 1 cos 



(161) 



where /o('2^) stands for the modified Bessel function. This expression provides us with the explicit form of the 
distribution of (scaled) partial delay times for the case of equivalent channels. Below we briefly analyze its most 
important features. 

First of all, we notice that the distribution above assumes the simplest form for the "critical" coupling T = 1 (i.e. 
g ~\) corresponding to the most strong overlap of individual resonances allowed for the few-channel scattering, see 
preceding sections. Under this condition one finds the following distribution of scaled partial delay times: 



-Af-2g-l/T 



(162) 



The powerlaw tail t^*^^^ at ^ 5 which is evident from the expression above for t; = 1 is actually a typical feature 
of the time delay distribution for any values of the parameters ga'iO. = 1, ...,M (see the discussion of Wigner-Smith 
time delay in the next subsection). For the equivalent channels we obviously have: 



(163) 



where Puig) stands for the Legendre polynomial |111]. For the case of non-equivalent channels the asymptotic 
behavior t~^~^ can be infered from Eq.(159) upon noticing that: Ak — J2b=i -^-b Ylc^bi^c — Rb)~^ = for fc = 
0, 1, Af — 2 and Ak ^0 for A; > M — 1, so that the integrand is proportional to 1/t^~^ . Combined with the factor 
in front of the integral it gives the desired behavior. 



In Fig. 3 we plotted the distribution Eq.(161) for the opposite case of weakly open systems (the regime of isolated 
resonances: T <C 1, i.e. g ^ 1). Under this condition when g ^ Tg the modified Bessel function can be replaced by 
its asymptotic expression valid for large arguments. Taking into account also that g — \/ g^ ~ I 



find that the distribution function Eq.(161) simplifies to the following form 



2s 



at 5 ^ 1, we 



d 



M 



(164) 



This expression is correctly normalized and plays the same role for the distribution of partial delay times as that play ed 
by the distribution in the issue of the resonance width distribution. It is necessary to mention that Eq.(164) is 
valid as long as l,rs <C g. At larger values of the behavior changes to that given by Eq.( |l63| ). It is interesting to 
note that in the parametrically large region (2g)^^ <C 1, <C <? one can neglect the exponential term in Eq.(164) and 
reduce this distribution to the following form: 



P.(r.) 



(2M-1)!! 1 
7m 



2*-fM! 



-3/2 



(165) 



This Ts behavior taking place irrespectively of the number of open channels is therefore the most typical feature 
of the partial delay times distribution for the regime of isolated resonances. The origin of such a behavior can be 
understood analysing t he g eneral expression for the Wigner-Smith delay times, Eq.(45) ( see a more detailed discussion 
after the equation Eq.( |l95| )). 

At Ts ~ (2(;)~^ the distribution shows a maximum at a value 'P(ts) ^ g and then is cut off exponentially at smaller 



Vr [js « (2g) 1) - ^^/f^, (25r,)(^+3/2) 



exp - 



2gTs 



(166) 



All these features are evident from Fig. 3. 
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Having at our disposal the exact distribution eq. ( |161| ) it is instructive to calculate the mean value and the variance 
of the partial delay times. One finds: 



2tt ^ 

WK'- 



2M{T-^ - 1) + 1 
M - 1 



(167) 



The first of these relations is quite well known |8^,^,^,^,Q . It shows that the mean delay time (r) is determined 
by the mean level spacing A of the closed system and the number M of open channels. On the other hand the 
ma gnitu de of delay time fluctuations measured by the relative variance of the partial delay times distribution, see 
eq.(167), is determined both by M and T. Generically, the fluctuations are the weaker the larger is the number of 
open channels M and the stronger is the coupling to continua: 1 — T ^ T. Let us also mention as an interesting 
feature the divergency of the time-delay variance at M = 1, which is a consequence of the powerlaw tail t~^^~'^ 
typical for the distribution Vr{Tg). 

Here it is appropriate to mention that recently two other groups of authors [^,118| addressed the question of delay 
time distribution by different approaches. Gopar, Mello and Biittiker ||89[] verified numerically an old conjecture by 
Wigner |119| concerning invariance of poles and residues of the A'— matrix under a certain set of transformations, 
provided there is only one perfectly open channel: M — 1] T — 1 case in our notations. When combined with the 
distribution of residues, this conjecture was shown to produce the time delay distribution for all three symmetry 



classes ( orth ogon al, unitary and symplectic), the result for unitary clas s just coinciding with Eq. (|l62| ) for M ~ 1. 
In the paper |118| Seba, Zyczkowsky and Zakrzewsky arrived at Eq.(162) for arbitrary Al after a set of shrewd, but 
uncontrolled manipulations with eigenvalues and eigenfunctions of AT— matrix. Actually, these authors suggested the 
following general expression for the scaled partial delay time distribution claimed to be valid for arbitrary value of 
the coupling constant 7 and E = 0: 



-7/r 



lFi [M,M + 1, 



7~7 



(168) 



where iFi[M, M + 1, z] is the confluent hypergeometric function. 

This expression is quite different from ours given in Eq.(161), the two formulas coinciding in the limit 7 = 1 only. 
However, one can check that the expression Eq.(168) fails to fulfil the following important condition: the mean delay 
time must be independent of the degree of coupling to continua, measured by 7. Instead, it should be determined by 
the mean level spacing A. This requirement is satisfied by our distribution Eq.(161), see Eq.(167) and is known for a 
long time ||8(]|,p2|,p3|,p5|j87|] . It follows from the basic formula (45): 



2N 

H 



-jy dT dup{E-uj-T)- 



r/2 



r2/4 



(169) 



where p{E; T) = ;^(X]^=i ^{-^ ^ En)5{T — r„)) is the density of the S— matrix poles. For few-channel case M N 
the typical scale of the width F is the mean level spacing A = 1/ {Np{E)), see Sec. HI, so that the Lorentzian factor 
in the integrand of Eq. (|l69|) ca n be replaced by 2Tr5{uj) when evaluating the integral. This gives (r) = 2ttNp{E)/M 
in full agreement with Eq.(167). At t he sa me time, for the particular case M = 1 one can find that the first moment 

in contradiction with the general 



corresponding to the distribution Eq.(168) is given by: ^{t)szz= 2 In 7/ (• 



(7-7 



discussio n ab ove. This failure rules out the distributon Eq.(16S) as the correct one and shows that the assumptions 
made in |118| are justified only as long as T = 1. 

The distribution Vw{w) of parametric derivative s of phase shifts Wa = '^^^ found in a very similar way. 
Proceeding in the same manner as in Eqs.(152-155) one obtains: 



Vwiw) = ^ (^^^^S {w - Wa)^ ^ J J dzduVx{z,u)S + Y 



2u 



(170) 



where 



rx(z,u) = j-A2^5{z~Za)5 \u 



(171) 



M lim XJCe n=Qx{zi ^ z- uX/2,Z2 = z + uX/2) 



Performing the limiting procedure AT — > in the same way as — > in Eqs. (153-156) we arrive at the expression: 
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I r dt ,2/2 92 M ^ 



loo (2^)1/2 ^ 2M^2 y_ ^2 11 1 „ ,^^^,-1 



which can be written as: 

rx{z,u,) = 1^ (2^)V2 (173) 

where = u/(ttVsc) and VEiz^Vs) is the joint probabiUty density of z and its derivative over the energy studied 
earher in this subsection. This fact means that a similar relation Eq.( |l73[) holds for the distribution function Vw{ws) 



of scaled parametric derivatives Ws = ir-^ — §^ and that of the scaled partial delay times Vt{ts)- 



00 



The same relation was obtained in the paper |llq| on a basis of some plausible assumptions concerning parametric 
derivatives of phase shifts. 

C. Parametric correlations of Wigner-Smith time delays 



The expression for the correlation function Eq.(14C) can be used to calculate the parametric correlations of Wigner- 
Smith time delays Tyj{E,X) defined as 

r^(£;,X) = ~^Aindet5(i?,X)^l:A^0^(i5,X). (175) 

a 

To show this, we should remember the relation Eq.(^^ between the total phase shift 6 = 0a, the exact density of 
states for the closed chaotic system vx{E) and the eigenvalues Za of the iiT— matrix. We see that 

2'kN , ^ , , , , 9 1 J- ^ , > , ^ 

T^{E,X)^^vx(E) + T,{E,X)- r,(i?,X) = -2— —^arctanz, (i?,X) (176) 

a— 1 

Introducing the correlation function: 

Cw{^,X) = (t,,{E - n/2,-X/2)T^{E + n/2,X/2)) 
we see that it consists of three essentially different contributions: 

Cwin,X) - Crr{^,X) + ^ iC,ri-^,-X) + C,r{^,X)) + {2tiN / Mf X) (177) 

where 

Crr{^,X) = {t,{E - n/2,~X/2)T,{E + n/2,X/2)) (178) 

C^ri^,X) = {iy{E -n/2;~X/2)T,{E + n/2,X/2)) (179) 

C^^{n, X) = {v{E - 0/2, -X/2)v{E + n/2, X/2)) (180) 
In what follows we are interested, as usual, in finding the "connected" part of all these correlation functions. This 



will be implicitly assumed below. The correlation function Crri^iX) can be easily related to that given by Eq.(14C) 
because of the relation: 



d f°° 

Tz{E,X) ^ ~2— / c?zarctan(z)pE.x(2:) 
OE 
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where pe.x{z) is the density of JC— matrix eigenvalues defined in the beginning of this section. As a result we have: 

92 Q2 



nOO nOO 

j — / dziarctanzi / arctan22/c('2^i, -2^2) 

J — 00 — 00 



(181) 



where we used the relation Eq.( |l26| ) between the correlation function IC^^ujizi, Z2) and that given by the Eq.( |l4C| ) 
and we used the scaled variable u — irfl/A. After such a rescaling it is obvious that the term containing the second 
derivative can be n eglected i n co mparison with the second one because of the large factor A~ . Substituting now 
the expression Eq.( |l40| ) into eq.(181) one can easily perform the integration over 2:1,2:2 by exploiting the presence of 
the second derivative d^/dzidz2 in the function /c(zi,Z2) (this allows to convert factors arctanzi.2 into (1 + zf 2)'^ 
by partial integrations) and noticing that all poles of the expression 



M 

n 



zi(a)z2(a) + ijaMiziia) ~ 2:2(0)) +7a"' 



21(0)22(0) +i7aAi(2i(a) - 22(a)) +7a 



lie in the upper half plane Im2i > with respect to the variable 21 and in the lower half plane Im2i < with respect 
to the variable 22- As a result, the integration can be performed trivially by closing the integration contour over 
21(22) in the lower (upper) half plane, correspondingly, and amounts to replacing 21 = — «; 2:2 = i in the integrand. 
This gives: 



Crr{cJ,X)^2(^j^j ReJ d\2 j dAiexp{zw(Ai 



-A2)-y(A2-A2)}/(Ai,A2) 



where 



M 



/(Ai,A2) = l- 



nl 1 + 7~aA2 - ilaE/2 



^ ( \+ia\2^llaE/2 \ 1 -i^ 

l}^\l+i.X,+^^.E|2)] n 



1 + 27,A2 + 7g 
1 + 27aAi + 72 



(182) 



(183) 



The first three terms are boundary contributions due to the partial 21, 22 integrations. 

Let us now show how to calculate the correlation function C^ri^, X) which is expressed to the leading order in 
iV > 1 as follows: 



d f°° 

Ci,t{^,X) ^ — J dza.vcta.nz{v_x/2{E -n/i)pE+n/2.x/2{z)) 



(184) 



To this end we represent the density vx{E) in terms of the corresponding resolvent as vx{E) — :^Iiii Tv[E 
iO+ — Hin)^^ so that 



1 



2'k'^MN 



Re ( Tr- 



1Cut{^,X;z) = {i^-x/2{E - n/2)pE+n/2,x/2{z)) 

1 nn 1 

-Tr- 



(185) 



E-n/2- H{-X/2) - ?:0+ 2 - K^Q., X/2) + ^0^ 



The calculation of the connected part of the averaged product of two traces of resolvents in the preceding equation (we 
denote this quantity henceforth as fc{z) omiting an explicit mentioning of th e pa rameters ^,X) goes along exactly 
the same lines as the calculation of the correlation function /c(2i, 22), see Eq. (|l25| ). Namely, one writes this function 
as: 



^ciz) 
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dJidJ2 



Z 



(2) 



M 



J-y{Jl, J2 



T>ei[E - iO+ - n/2 - H{-X/2) + Ji]Dct[S + n/2 - H^{X/2; zf^)] 
Det[E - iO+ - n/2 - H{-X/2)]Bet[E + VL/2 - H^{X/2] Z 



where we introduced the notations: Z 



(2) 



^0^ 



J2 and H^{X;Z 



(2)^ 



Hi; 



(2) ^1 
X 



H 



(1) 



(186) 



-^WW'f. The 



^ J 

generating function J-^(Ji, J2) is expressed in a standard way as a Gaussian superintegral. Finally, the function f^{z) 
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is reduced (after exactly the same manipulations as before, see Eqs.( 12^ 134)) to the following representation in terms 
of the nonlinear cr-model: 



r (11) f ^ \ *^ 

fciz) = X / - ^yQ\ [ n Sdet-iB,(a) ) ^ {B-\a)) 



(22) 



(187) 



exp 



-^^StrQA + ^StrQAQA I 
2 8 1 



where the supermatrix By{a) is given by 



B,{a) = 



1 + A 1 - A 



(^z(a)/4 + 7a(3 



and we used the conventions of Appendix B for matrix elements as well as the notations: z{a) = z — 7ai5/2; 7q = 
T^i^scla- The exp licit expressions for all matrix elements are given in Appendix B. Substituting them into the super- 
integral Eq.(187) we find again, that the term containing no Grassmannians at all gives the contribution 



M 



N{E/2 + iTTv)^- 



1 



^ z{a) + ila 

which is exactly the "disconnected part" of the corresponding correlation function. The connected part is given by: 



A2 



^— ^ exp{iw(Ai - A2) - y (A^ - A^)}— I I ^ '— 



z ^\ z{a) + i7Ai 

a— 1 ^ ^ 



(188) 



Substituting this expression into Eq.(184) one can again trivially perform the integration over the variable z. As a 
result one obtains: 



C^^(cj,X) = - ^J^jj^ (^) Re y" d\2 dAiexp{iw(Ai 



M 



A2)-y(A?-A^)}x 



1 + 7aA2 + ^7a£^/2 
+ faAi + i-iaE/2 



(189) 



Finally, the parametric correlation function (7^^,(17, X) of the densities of states for a closed chaotic system with broken 
time-reversal symmetry was found some time ago by Simons and Altshuler [E8|: 



„{a,X) = — (-j Re j ^dX2 dAi exp{ic^(Ai 



A2)-y(A?-A2)} 



(190) 



Summing up all the contributions we find the desired expression for the parametric correlation function of scaled 
Wigner-Smith time delays fg 



A , 



2M2 



dX 



Cw{uj,x) = {5fw{E -n/2,-X/2)Sfw{E + n/2,X/2)) = 

QQ. 2 ^ r -1 \ 1 

dAiCosKAi-A)]exp{-^(A?-A^)}n T^T^ 

2 „V;U + Aiga 



(191) 



where we used the parameter ga = 2T^^ — 1 introduced earlier. 

It is interesting to mention that there exists an alternative way to derive the pair correlation function of Wigner- 
Smith time delays given in Eq.(191). The starting point in that case is Eq.(^). Then the calculation of the correlation 
of fluctuations of Wigner-Smithtime delay 5tw{E, X) — t^j ~ {t.^) amounts to evaluating the average product of the 
resolvents of the non-Hermitian effective Hamiltonians Ti, ± iirWW^ . This can be done by exactly the same method as 
we use elsewhere in the present paper. For the case of chaotic systems with preserved TRS and no external parameter 
X such a calculation was done earlier in [ pTf . 

L et us analyse the correlation function Cwix^uj) in more detail. For this purpose we find it convenient to rewrite 
Eq.(|l9l[) in a slightly different form: 
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where 



/•oo 

R^^''^\uj,x) = / dt I COS {uti) exp 
Jo 



-xHi 



xHi 
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TT — — 



/2 



(192) 



(193) 



dt2 cos {ujt2) exp 



Y[il-TJ,/2) 



(194) 



and the functions R^^j'^\uj,x 



P 



1,2 are obtained from the expressions for R^^j'^\lli,x); p — 1,2 by replacing 



cosujtp by sincjip. 

First of all, the correlation function Eq.( |192| ) taken at uj = x — gives the variance of the Wigner-Smith time 
delay distribution. In principle, the corresponding integration can be performed for an arbitrary set of transmission 
coefficients Ta- The resulting expressions turn out to be quite cumbersome. They simplify in the case of all equivalent 
channels T,, — T: a = 1, M when we find: 



T,:, - r,. 



T2(Af2 _ 1) 



(195) 



0) as those following from 



Thi s exp ression shows the same qualitative features ( divergencies at M = 1 or T 

Eq.deJ. 

For chaotic systems with only one open channel the Wigner-Smith tim e de lay fw just coincides with the partial 
delay time and the corresponding distribution is given by that in Eq.( |l6l[ ). Unfortunately, our methods give us 
no possibility to find explicitly the distribution 'Pw{t'w) of Wigner-Smith time delay for an arbitrary number of 
open channels M > 1. It is natural to put forward a conjecture, that the divergence of the variance of Wigner- 
Smith time delay as long as M — > 1 indicates that a (unknown) distribution Vw{tw) possess the same powerlaw tail 
Vwi^w) oc T^~^^~^ at large fw as that typical for the distribution of partial phase shift times. Another argument 
supporting this conjecture comes from the general formula Eq.(45). Taking the value E at random it is evident that 
anomalously large time delay Tw{E) ~ F"^ corresponds to the event when E happens to be sufficiently close (at the 
distance 5E ^ F„) to a position En of an anomalously narrow resonance F„ <C A. The probability of such an event can 
be estimated as Pr oc (r/A)p(r/A <C 1) o c (F/A)^^, where we used the small width asymptotic p(y <C 1) oc of 



the resonance widths distribution Eq.( 102 ). Then the asymptotic tail of the probability distribution of the time delay 



can be estimated as 7'(tu,) cx J dT6{Tiu — F ^)Pr oc r^, ^ in agreement with our conjecture. Here it is appropriate to 



mention that the same asymptotic behavior is typical for the staying probability function p{t), see Eq.(p6h. The long- 
time asymptotic for p{t) was found for the systems with preserved time-reversal symmetry in the papers | |25| , |8^ , p6[ . 
It is trivial to adjust the corresponding argumentation to the present case and to recover the behavior. 

The expression Eq. ( p5[ ) allows one also to show that for weakly opened systems the distribution of the scaled delay 
times should demonstrate the universal behavior V{t) cx r"'^/^ in the parametrically large domain <^ t <^ g , ci. 



Eq.(165). Indeed, for 5^1 the resonances do not overlap: F„ <C A, and their widths F„ follow the x^distribution. It 
is therefore clear, that for any particular value of the energy E the sum in Eq.(^5|) is dominated by a single resonance 
whose position i?„ is the closest to E, that is 7\„ w Tty'^+v? ' '^^ere t/„ = ttF/A and u„ = ^{E — En)- We can 
estimate the value of the contribution coming from all neglected terms assuming that all resonances have the same 
widths (F) and are spaced equally with the mean spacing A. This immediately gives the correction to be of the order 
of 5fw oc (F)/MA ^ g~^, where we used the formula Eq.(llC) in the limit T ^ 1. We conclude that the distribution 
of the scaled time delay is correctly reproduced by the distribution of the" closest resonance term" as long as we are 
interested in the region ^ 9^^- Assuming that the variable w„ is uniformly distributed in the interval [— 7r,7r] we 
find: 



Vif^)= / dyr^2{y) / du6 



y 



dy 



y 



1/2 



(2- 



(y) 



(196) 



where ym{fw) — min ; tt^ ^'^^^ . Taking into account that the distribution V^^iy) is cut exponentially at 
y > 1/g we can safely neglect the term Mry <C 1 and set the upper limit y„i = oo as long as <^ g. This 
immmediately results in the anticipated f^^^'^ behavior. On the other hand, we can put V~^2{%j) oc g^'^y^^^^ and 
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Um = 2/{Mfw) in the domain of extremely large time delays f^, ^ g, which results in the g^f^'"^^'^^'' tail in full 
agreement with the general discussion presented above. 

The behavior of the delay time distribution in the domain of extremely small time delays f^^ < is determined by 
contribution of many resonant terms in the expression Eq.(Ma). However, one can argue that th e dis tribution should 



be exponentially cut: ^'(t^ oc exp— ^const^rg) as is indeed seen from the expression Eq.(166). This behavior 
is a typical one for a sum of random variables of the form ^„ y„/u~^, with (j/„) ^ g^^ (the so-called stable Levy 
distribution, see similar arguments in |]66| ). 



The correlation function Eci.(192) acquires quite a simple form for the case of many weakly open channels: Ta <C 
I, a = 1,2,.., A/ but r = Y^a^a > 1. Then we can put effectively: - t2Ta/2) w exp-ri2/2; JlaCl + 

T 

corresponding integrals can be calculated exactly giving: 



tiTa/2) ^ Ri exp— rti/2 and also neglect the terms ± — in the exponents of the integrands in Eqs.( |l93|Jl94| ). The 



(r2^-cj2)[l-e-2rxcos2cj] +2rxwe-2rxsin2w 
Cw{^,x)^ (^2^p2^)2 (197) 

where Tx = T + x'^ . Let us note that for x = this expression is actually valid for arbitrary F. For arbitrary value of 
X the condition of validity is F^ 3> 1. Neglecting the exponentially small terms we arrive at the simple expression: 

^^(^'^)= (3^^^ (198) 

As will be shown in the next subsection, this expression is nothing other but the semiclassical formula for parametric 
time delay correlations in systems with brok en ti me-reversal invariance. For the case of preserved time-reversal 



symmetry and no external parameters the Eq.(198) was derived in l8 



D. Semiclassical theory for parametric correlations of time delays. 



A general semiclassical expression for the Wigner-Smith time delay in terms of a periodic orbit expansion has been 
given by Balian and Bloch [g7|. It is formally identical to Gutzwiller's trace formula. The corresponding expression 
for the pair correlation function of time delay (without taking into acco unt a parametric dependence) for chaotic 
scattering was derived by Eckhardt |2^. In parallel. Berry and Keating |120|| developed a method allowing to take 
parametric correlations into account for the case of a closed chaotic system pierced by a magnetic flux serving to 
break down the time-reversal symmmetry. Below we show briefly how to combine both approaches to arrive at the 
semiclassic expre ssion for the parametric correlation function of time delay in that case; see also related discussion 
in the papers ||lC , 121 , 122 1 . The semiclassical periodic orbit expansion for the "fluctuating part" of a time delay of a 
quantum particle with an energy E + f2/2; Q <ti E moving in a systempierced by a magnetic flux line with flux 
(measured in units of flux quanta 0o = 27rc/e) is: 



6T^,{E + n/2,(j)) = ^v4j- 



(199) 



where the summation goes over all periodic orbits with period Tj = dSj/dE, with Sj{E) being the corresponding 



action, A, 



being the amplitude and /Uj , Mj being the Maslov phase and stability matrix corresponding 



2iT--\/det(A/j-i) 

to the given periodic orbit. The winding number wj counts the number of times the orbit winds around the flux line. 

Thus, for the parametric correlation function one finds: 



Cw{^,X) = {Sfw{E -n/2,(l>- X/2)Sfw{E + n/2,<j) + X/2)) = 

mX{wj +Wk)\ ) 



(200) 



^ I A,Ak I exp ^ J [S,{E) ~ Sk{E)] + ^(T, + T,) -f 2^i{w 



2h' 



where the averaging goes over the energy spectrum. According to standard argumentation [22 120| one can restrict 
oneself to the so-called "diagonal approximation" taking into account only contributions with coinciding indices j = k: 



C#"^(f7,X)= 5]|A,pexp 



■ n 

i—T-j + 2T:iXw-j 

h 



(201) 
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The next important step uses the fa ct that w inding numbers for orbi ts in any narrow window of periods are essentially 
irregular and Gaussian distributed 123 12Cf| , see also discussion in 1 121, 122 1: 



where the variance (j{Tj) increases linearly with period: (j{T) — (3T. The constant [3 is system dependent; for a 

particle with mass m moving in a billiard of the area A it is proportional to {2E /mA)^^'^ | 122| , 121 |. Taking the 
discrete nature of the winding numbers into account one can write: 



exp 2'KiXwj 



Eoo 
k=- 



,-(2x)V(T_,)(X-fe)^ 



Eoo 
fc=-c 



,-(2Tr)2cr(T,)fc2 



(202) 



where we used X <C 1 (i.e. change of the magnetic flux is much smaller than 0o) and neglected exponentially small 
terms O (exp — 47r^cr(r)) . Substituting this average in the correlation function of time delay it is convenient to consider 

its Fourier transform C{t,X) = J dQe-'^^/'^C^'''^ (n, X). We have: 



Cw{t, X) 



E 



I A, p S{t - T,)e 



(203) 



where we denoted x = 2ti(3^''^X. 

For closed chaotic systems the sum in the preceding equation is known to be proportinal to the time t (this is the 
famous Hannay-Ozorio de Almeida sum rule |124|). For open chaotic systems Eckhardt gave some arguments in 
favor of replacing this sum by te'^"'*, with F^/ being the classical escape rate from the chaotic region. Using this fact 
we see that (under the assumptions we made) the semiclassical expression for the Fourier transform of the correlation 
function of time delays is given by: 



Cwit,X) ^te 



(204) 



Afte r Fo urier-transforming this expression back we see that the result turns out to be identical to that given in 
Eq.( |l98| ) upon identification V^i — T; i — > a;. 

Another interesting point to be mentioned is that the form of the time delay correlation function given in the 
eq.( |197] ) (which contains Eq.(198) as a limiting case) was obtained by Shushin and Wardlaw |125| in the model of 
chaotic scattering on a leaky surface of constant negative curvature. At the first glance such a correspondence is 
quite a su rpris ing fact since the model considered in [125| corresponds formally to one-channel scattering, but the 
result Eq.(197) was derived under the assumption of many weak chanels. In order to understand that fact one should 
remember that the model considered in |125{ possesses quite a peculiar property: all its resonance poles turned out 
to have exactly the same widths. It is at variance with the known form of the resonance widths distribution for one- 
channel scattering in a generic chaotic system, see eq.( |102| ), where resonance widths fluctuate strongly. At the same 
time, if we consider the limiting case of many weak channels: M ^ 1,5 ^ 1 and M/g = T^f fixed , the distribution 
of r eson ance widths tends to the delta-functional one p{T) — 5{T — T^f). This fact can be easily infered from the 
Eq.( |l08[ ). We see that effectively it is just the limiting case of many weak channels that corresponds to non- fluctuating 
resonance widths. Under these conditions the correlations of the time delays are determined by the statistics of the 
positions of resonances. For the model of scattering on a leaky surface of negative curvature the positions of resonances 
are given by the zeroes of Riemann zeta-function on the so- calle d critical line in the complex plane. Accordi ng t o 
the celebrated Montgomery conjecture (verified numerically |126| and supported by sound analytical results |127|) 
statistical properties of these zeroes are identical to those of eigenvalues of large random GUE matrices. All these 
facts taken into account it is no more a surprise that the correlations of time delays for both models coincide in the 
considered region of parameters. 



VI. SUMMARY AND CONCLUSIONS 



In the present paper we analyzed in much detail the universal features of statistics of resonances, phase shifts 
and delay times for a generic open chaotic quantum system with broken time-reversal invariance. This was achieved 
by replacing the Hamiltonian of the chaotic region by a large Random Matrix taken from the Gaussian Unitary 
Ensemble. Employing the well-developed method of mapping the problem to the so-called supersymmetric nonlinear 
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cr— model we succ eede d in d e rivin g expl ici t an alytic al ex pressions for various distributions and correlations functions 
, see Eqs.(|97|-p9|), (|l02| ),(140),( |l59| ),(161), ( |l74|) and (191), characterizing the above mentioned quantities for arbitrary 
finite number of open channels and arbitrary strength of coupling to continua. 

The best candidates for checking the validity of t he e xpressions obtained are realis tic models of mesoscopic ballistic 
devices subject to applied magnetic field. Closed |128[ as well as open |^,122,12S| systems of this kind were inten- 



sively investigated recently and the statistics of S-matrix elements and related quantities was available among other 
char acteristics. Very recently, the issue of dwell times inside the chaotic region attracted some research interest as well 
1 130 1 . All these facts allow us to expect that our results can be verified independently in the numerical experiments. 
It is interesting to mention that recently another type of chaotic systems with broken time-reversal invariance became 
available experimentally |l^ . The authors used microwave resonators of billiard shape with a " handle" which allows 
only unidirectional propagation of radiation in it thus breaking the symmetry between the wave and its time reversal 
counterpart. 

It is important to mention that our results are also of potential experimental relevance. Indeed, the issue of time- 
delay fluctuations turns out to be intimately related to the statistical properties of mesoscopic capacitors p8| , p9[ . 
Considering the case of a mesoscopic cavity coupled by a M-channel lead to one electronic reservoir and capacitively 
to another reservoir Gopar, Biittiker and Mello |^ suggested the following expression for the low frequency AC 
admittance of such a structure: 



(205) 



where Ce stands for a geometric capacitance relating the charge Q on the plate to the voltage U across the capacitor, 
Tw stands for the dimensionless Wigner time delay, and rj — j^§fj-^, with A being the mean level spacing for the 
cavity. 

For macroscopic cavities — > and the dimensionless capacitance a is equal to unity resulting in the classical 
expression for the capacitive response: {cu) = —iuuCe- In contrast, for small enough cavities rj has to be taken into 
account and the fluctuating delay time t^, results in a fluctuating admittance. 

As it was discussed above, for one open channel M — 1 th e di stribution of Wigner time delay is identical to the 
distribution of partial delay times Vt{t) and is given by Eq.(161). This fact immediately allows one to write down 
the distribution of the dimensionless capacitance a as: 



V 



(1-ay 



rja 



(206) 



For the perfect coupling case T = 1 the corresponding distribution was analyzed in ||89| . In the opposite limiting case 
of weak coupling T ^ 1 the universal t~'^I'^ time delay distribution, see Eq.(165) results in the following expression: 



-1/2 



-3/2(1 -a)-i/2 



(207) 



as long as a 3> Tjr] and \ — a rjT. As it follows from our previous discussion after Eq.( |l95| ), this form of the 
distribution should be valid for arbitrary number of weakly open channels. 

Actually, our knowledge of the general expression for Wigner time delay variance , see Eq. ( |195| ) , provides us with 
the possibility to determine the variance of the low-frequency admittance G^ [oj) in the limit of many open channels. 

Indeed, in the limit Af ^ 1,T ~ 1 our expression just says that the variance of time delay is of the order 
l/A/2 ^ 1 as compared with the squared mean value (r„)2. Thus, we can represent the fluctuating time delay in 
a form t„ = (r) -I- St ^,, w here typical scale of the fluctuating part is of the order of St^ ~ {tw)/M. Substituting 
this expression to Eq. (|205| ) and expanding with respect to the fluctuating part 6t^ one obtains to the first nontrivial 
order: 



G' {io) ^ {G\uj)) 1^ + 57 



(tw) (»7+ (r^)) 



and immediately extracts the variance of the admittance : 

{{G'iu;)f)-{G'iu;)y 



' RC 



T2Af2 (r,„)2 



where = {t^)~^ + is the so-called RC-time and we substituted the expression Eq. (|l95| ) for the time delay 
variance taking into account that A/ 3> 1. Such an expression for the particular case T = 1 was very recently derived 
by Brouwer and Biittiker by a different method M . 
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The majority of the numerical data concerning various statistical properties of the sc attering matrix for open 
chaotic systems corresponds to the case of preserved time-reversal invariance see e.g. 72 131,132]. This case is not 



only simplier from numerical point of view ( opposite to the situation with the analytical calculations), but also the 



most relevant experimentally. As a consequence, numerical studies on resonance width statistics [|16| , |17 52 
as well as on properties of scattering phase shifts and their derivatives [ f75| , [74|j76| were restricted to the systems of 
that symmetry class. It is necessary to note that some analytical results for poles and time delays for time-reversal 
invariant scattering are already available in the literature for some time. In particular, for only one open channel the 
joint probability distribution of all N complex resonance poles is known |]53| ( however, not the density of these poles 
in complex plane) as well as the distribution of Wigner time delay for perfect coupling to continuum . Essential 
progress was achieved by Lehmann et al. |Q who calculated the correlation function of time delays for two different 
values of energy and any number of open channels. Actually, the calculation similar to that done in the present paper 
can be successfully carried out for the whole crossover region between the orthogonal and unitary symmetry classes. 
The results will be published elsewhere 

Let us also mention that recent numerical results jOSt show that the resonance widths distribution derived in the 
present paper can be applied for the systems with preserved time reversal invariance quite satisfactorily after replacing 
the number of channels AI by Af/2. This fact is not so surprising, taking into account that such feature as the powerlaw 
tail 1/y^ of that distribution is actually a generic property following from the chaotic classical dynamics only, see the 
discussion in Sec. III. 

For a majority of models in atomic and molecular physics parameters of all the resonances can be determined even 
without expensive calculations of 5— matrix elements. The most effective method is the so-called complex scaling (or 
complex rotation) method |^l[] successfully used for the systems exhibiting chaotic behavior JT6| , pT|j5^ ] . It is interesting 
to mention that a crossover from isolated to overlapping resonance regime was detected recently for the dissociation 
reaction HO2 — > i? -t- O2 in one open channel case |6^. One can hope that applying the complex rotation method to 
this sytem one could extract the widths of resonances with sufficient accuracy and to observe a transition from the 
distribution towards that with the l/y^ tail. We would like to point out that the whole 5- matrix as a function 
of energy of incoming waves was measured in real experiments |9|,[l0|, and even used to calculate the average time 
delay ||ict| . In principle, the positions of resonances in the complex plane can be extracted if one knows S{E) with 
sufficient accuracy. For example, one can use that fact that for any number of open channels the determinant detS{E) 
as a function of energy has its singularities (which are just resonance poles £'„ — «r„/2) only in the lower half plane 
Im£ < 0. As a result it can be written as: 

where 5 is the phase of potential scattering irrelevant for our discussion. Provided the values of detS{E) for real E 
are known, one can restore the determinant of S"— matrix in the upper half-plane lm.E > by the relation: 

detS{E + il) = - / -^^ -jdet5(£;) (209) 



It is easy to see that the two relations ( p08D and (20£) allow to determine all the resonance parameters i?„, r„ from 
zeroes of the S'— matrix determinant in the upper half-plane Im£ > 0. Of course, the practical implementation of this 
procedure requires highly accurate data for S{E) which is not the case in the mentioned experiments due to noise and 
damping in resonator walls. However, one can hope that the progress in the experimental set-up could make such a 
measurement feasible in future. 

Finally, as an interesting perspective for future research we would like to mention the issue of S"— matrix statistics 
for systems exhibiting the Anderson localization phenomenon. This issue attracts research attention for some period 



[133| and increasing amount of numerical results are already available |75,11C ,7^,134| requiring a systematic analytical 



insight into the problem. 
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APPENDIX A: REGULARIZATION OF THE EIGENVALUE DENSITY FOR NONNORMAL RANDOM 

MATRICES 

Let us consider a random, but fixed non-Hermitian N x N random matrix TC, for which we only assume that 
generically its complex eigenvalues are non-degenerate. Taking the second derivative of the potential Eq. (p9[) ( apar t 



from the factor with respect to the energy £ and its complex conjugate £*, we obtain Poisson's equation [[135 

1 1 

" dsds^'^ = (H - sViH -£) + K^ in -E){n- f )t + k^= ^^^^ 

with a density which is always positive, because the operators appearing are both positive. We will show below 
that Pk goes to a sum of two-dimensional (5— functions in the complex energy plane: 

N 

^^^p,{E,Y)^Y.^^ (A2) 

where £j are the eigenvalues of Ti. The weight for each (5— function is one. Here the integral over p{E, Y)dEdY is 
normalized to N: J p^{E,Y)dEdY = N, which can be kept finite. This is true for any k > and can simply be 
shown by Stokes theorem. 

Now let us consider the Hermitian eigenvalue problem 

[in - £y {n-£) + «2] = \,^, (A3) 

The eigenvalues we again assume to be generically non-degenerate, and the ipi form a complete orthonormalised 
set. It follows: 

[{n -£){n~ £)^ + k"] [n - £Mh ^K{n- (A4) 



so that (pi — {n ~ £)il}il — is a normalised eigenvector of a second Hermitian eigenvalue problem for the 
operator [(H — £){n — + k^] . Such an eigenvector corresponds to the same eigenvalue Xi, provided Xi ^ k^, and 
£ is not an eigenvalue of n. If {n — £)ipo — 0, we nevertheless can find a normalised eigenvector ipo, orthogonal to all 
i 0, with {n — £){n — £)^ 4>o = 0, i.e. 0o is an eigenvector of with the eigenvalue £* . 
Now we may expand in terms of these eigenfunctions: 

p.(i?,r) = ^^l|(^„0fe)p^ (A5) 

ik 

from which one sees explicitly that is positive. Here (i/;^, (j)k) stands for the complex scalar product in Hilbert space. 

If £ is not an eigenvalue of n (and £* is not one of n^), then goes to zero proportionally to for k 0. 
Therefore the weight to the normalization of comes only from the neighborhoods of the eigenvalues £j of 7i. If £q 
is exactly an eigenvalue of Ti, then we know the lowest eigenvalue Aq of {n — £q)^{1-1 — £q) + k'^ { which is obviously 
Ao = K^) and all other eigenvalues are higher by amounts independent of k. That means: 

p,(i?,r)« ^ I (V'0,'^0) P (A6) 

for £ = £0 and k — > 0, which diverges as it must for a (5— function at an eigenvalue £o- 

In order to see how p^ varies with energy near an eigenvalue £0, we may set £ — £q + S£ and calculate Aq by 
perturbation theory. Only the second order perturbation contributes and the surprisingly simple result is: 

Xo^K^+\6E\^\{^o,^o)\^ (A7) 
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This means that in the neighborhood of an eigenvalue £q the function has the form: 



P^^^^^>-- in^+\6E\^\{4,o,i.o) ^ ' 

We will not consider rare cases, in which the vector V'o is occasionally orthogonal to (/)o. Then we see that Pk.{E, Y) 
has in the limit k ^ the form of a two-dimensional 6— function with weight one. Its width goes to zero like 
k/ I ((/)o,V'o) I ■ T his is valid in the neighborhood of one isolated eigenvalue, which is however arbitrary. This proves 
that the Eg. ( |A2D is indeed correct. 

APPENDIX B: THE PARAMETRIZATION OF THE MATRICES Q 

The supermatrices Q belonging to the graded coset space f7(l, l/2)/J7(l, 1) x U{1,1) can be parametrized as 
Q = XJ-^MU where 

Ml I e*^^ \ ,r -f \f^i\ 



where a, a*, (3, (3* are Grassmann variables, Ai € (1, oo); A2 G (—1, 1); 4>i, (1)2 G (0, 27r) and Ai 2 are related to | /ii 2 
via Af- I ^1 |2= 1; A2+|^2p=l. 

It is convenient to have also the explicit expressions for the matrix elements ofQ^f^^^ /^^^)■ We have 



Q21 Q22 

i/Ai \ / -i[Ai -a*a(Ai - A2)] ia*(Ai-A2) , , . 

[0 X2j''-[ ia{X,-X2) -t[X2-a*a{Xi~\2)] I 

- I A2 ) " = i /3(Ai - A2) ^ [X2 + /3*/3(Ai - A2)] ) ^^^^ 



„ _i / ui 
Qi2=^i ^^,*2 

^1(1 - a*a/2)(l + 13* (3/2) - a*Pfi*2 -i(3*{l - a*a/2)^i + - /9*/3/2)^* 

-a(l + (3*P/2)ni ~ /3(1 + a*a/2)^^ + i{l + a*a/2)(l - [3*(3/2)^l*2 



Q21 = V \ \{ ■ ] u = 

filil - a*a/2)(l + /3*/3/2) + a/3>2 -/?*(! + a*a/2)fi2 -a*{l+ /3*/3/2)^*i 
- (3*(3/2)fi2 - i(3{l - a*a/2)fil ~ia*(3fil + i{l + a*a/2){l - (3*P/2)fi2 



(B4) 



(B5) 



where we introduced the notations: fii =\ /ii | e'*^^; /i2 =| fJ-2 \ e"^^. The expressions above are frequently referred 
to as the "Efetov paramctrization" for the matrices Q. 

We denote the corresponding measure as dQ. Straightforward, but lengthy calculation gives: 

= 7T — ^ ■ da* df3* dadl3 (B6) 

(Ai - A2)^ {2ttY 

In the rest of this Appendix we present the explicit expressions for supertraces, superdeterminants, and matrix 



elements, entering different expressions in the main text, see the Eqs. (p6, 135, 187) 
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1. For resonance widths calculation. 

We have: 



StrQA = (q(;i) - Qif - (Q(.y) - Q^f) = -2z(Ai - A2) 



StrQ(A - KbL) = Q''ff - Q^^f ^ 2i\2 - - \2){a*a - P* (3) 



StrQE.=z(Qr-Qif)-(Q(-)+Qf/ 



lif^i - Ait - - /^2) - 1 7, (a^i - + + A*2) - «a*^(M2 - Mi)- 



ia(3*{fii + ^2) - -Q;*a/3*/3(/ii - - fJ-2 ~ ^2) 



The superdeterminant Sdet 



can be easily evaluated because the supermatrix [/, 



Eq.(Bl), commutes with A and therefore can be omitted under the sign of the superdeterminant: 



Sdet' 



Sdet" 



(B7) 



The supermatrix A 



I + i\^aEh. + inVsclaQ^ 



is however block -diagonal in the fermion-boson arrangement: 



Af,f = Afh ~ 0, see Eq.(Bl), and therefore Sdet A — DetA/^/Detylbb. Trivial calculation gives: 



Sdet" 



l + 2TTV,c{E)-faX2+ll 
l + 2TTV,c[E)-fa\l+ll 



which is reduced to the form used in the text of the paper, see 

e.g. Eq.(p6|) upon introducing the transmission coefficients Tq, Eq.(|49|) and the parameters ga — 2/Ta— 1. 

2. For scattering phase shifts statistics. 

First of all, using UK — AU one has 

Str(QAgA) = 2Str(A/i\ - M12M21) = -4(A? - A^) 
The main object entering the calculation of eigenphases correlation function is the supermatrix ^^a^ , where: 

Ba = i^i(/4 + A) + iz2(/4 - A) - 7, (f A - ttuscQ) = U-^b{a)U 
b{a) = diag (^zi(a)/2, Z2(a)/2) + %M 



(B8) 



where we used the notations: Zp{a) — Zp — ja-^', P — ^,'2 and — T^i^scja, the supermatrix M being defined in 



Eq.(B2). One can invert ba easily noticing that in the boson-fermion arrangement this matrix is block-diagonal: 
b{a) = diag ( hb (a) , iff (a) ) , where 



htia) = ( 7^'^"^^ , f ^^-^^ , 1 and bff{a) ^ ( "^^^^ ^^Ta 

' " -2(a) + l^laM I ■'^ 



Ml 7a Z2( 



SO that 5 (a) = diag b (a) , b (a) 

J bb 



^J■2la Z2{a)+i^a>^2 



ff 



where 
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5(a) 



bb Vbia) \ -Mi7a zi{a)-i%Xi 



(B9) 



b{a) 



1 / Z2(a)+«7aA2 -^^J■2la 



ff 'Df{a) \ ~ifJ^2la zi{a) ~ i%\2 



and we introduced notations: 



X'fc(a) = detStb(a) = Zi{a)z2{a) + ija {zi{a) - 23(0)) Ai +7^ 
I?/(a) = det6//(a) = zi(a)z2(a) + Z7a (^i(a) - 22(0)) A2 + 7^ 



(BIO) 



(Bll) 



so that Sdet [b-^) = Sdet (6-^(a)) 



'Dbia) ' 



Rearranging the supcrmatrix b ^(a) in advanced-retarded order we can find easily the supermatrix 



U b {a)U. ActuaUy, we need only its elements in the fermion-fermion block (see Eq.(135)) 



Str (s-^Ci) = - (B-i)"^ ; Str (3-^02) = (B-^) 



ff 



and Str 



We find, correspondingly: 



nil _ Z2{a) + ijaX2 



ff 



a a 



Z2{a) + i%Xi Z2{a) + ijaX2 



Vb{a) 



(B12) 



/ 1n22 _ Zi{a)-i^aX2 

>ff^ Vf{a) 



zi (a) - z7aAi zi (a) - A2 



2?fc(a) 



(B13) 



(B-^' 



ff 



Vb{a) T^f{a) 



(B14) 



(B15) 



3. For time-delay correlations. 



The main new object here is the supermatrix (see Eq.( 187 ) 



BM 



1+A 1- A 



z{a)Ii + 7q(5 



U 



h 

7aM2i z(a)/2+i7aMi 



(BI6) 



where we again used that matrices U and A commute. The matrix B^{a) is simple to invert. Performing the 
calculation we find: 



(Buia)-' 



1 



1 



'// 2(a) + i7aA2 \z{a) + i^aXi z{a) + i^aX2 

and the corresponding superdeterminant is given by: 

z{a) + i7qA2 



Sdet" (a) = 



z{a) + i7aAi 



(BIT) 



(B18) 
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APPENDIX: FIGURE CAPTION 



• Fig 1. A generic model for chaotic scattering : an irregular shaped cavity attached to the perfect lead. The 
Hamiltonian associated with the cavity region is simulated by a random matrix Hin. 

• Fig. 2 The distribution of scaled resonance widths p{y) for M — 1 (solid) M = 2 (dash-dotted) and M = 3 
(dotted line) equivalent open channles. The effective coupling is maximal: g = 1. As the result, the distributions 
demonstrate M/ (2j/^) asymptotic behavior at large y. 

• Fig.3 The distribution of scaled partial delay times ■P(ts) for M = 1 (solid) M = 2 (dash-dotted) and M = 3 
(dotted line) equivalent open channles. The effective coupling g = 10 corresponds to weakly open systems. 
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